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Preface 



Lately, after ~ 50 years, there seems to be a convergence in the languages used to 
describe Quantum Field Theory and String Theory, such that it seems possible to 
relate objects from these two perspectives. Therefore, a deeper investigation of the 
properties and features of QFT is a reasonable thing to do: non-perturbativc effects, 
dualities, emergent properties, non-commutative structures, etc. This particular 
line of research uses Symmetry Breaking in order to probe a few of the different 
topics mentioned above, i.e.. Symmetry Breaking is used as an "underlying principle", 
bringing different features of QFT to the foreground. However, the understanding 
of Symmetry Breaking that is used here is quite different from what is done in the 
mainstream: Symmetry Breaking is understood as the solution set of a given QFT, its 
vacuum manifold, or, more modernly, its Moduli Space. Distinct solutions correspond 
to different sectors, phases, of the theory, which are nothing but distinct foliations 
of the vacuum manifold, or points in the Moduli Space (for all possible values of 
the parameters of the theory). Under this framework, three different problems will 
be attacked: "Mollifying QFT" , "Topological Transitions and Geometric Langlands 
Duality" and "Three-dimensional Gravity and its Phase Transitions" . The first makes 
use of the Modufi Space of the theory in order to construct an appropriate moUification 
of it, rendering it viable to simulate a QFT in Lorentzian spaces, tackling the "sign 
problem" hcads-on. The connections with Lee- Yang zeros and Stokes Phenomena will 
be made clear. The second will show that each different phase has its own topology 
which can be used as Superselection Rule; moreover, the Euler Characteristic of each 
phase gives it quantization condition. The mechanism via which several dualities 
work will also be elucidated. The last one will generalize a 0-dimensional QFT, via 
dimensional construction through its D-Module, and conjecture several connections 
between the Lie- algebra- valued extension of the Airy function and the recent Partition 
Function found for three-dimensional gravity with a negative cosmological constant. 
These three problems, put together, should exhibit a solid and robust framework for 
treating QFT under this new paradigm. 
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Chapter A 

Motivations and Introductory Remarks 



The paradigm of research in Physics and Mathematics around the 1960's was that 
Physics dealt with local structures (QFT, "Dual Models", etc) while Mathematics 
dealt with global tools (Morse Theory, Cohomologies, Index Theorems, etc). 

Nowadays, this situation is somewhat reversed, with physicists worrying about ques- 
tions of Yang-Mills over Riemann surfaces, Higgs Bundles, Topological Quantization, 
etc; while mathematicians have gotten their hands full with structures inspired by 
QFT and String Theory, e.g., D- Modules. 

Thus, it seems to be that there is a convergence of these pictures into a single 
language: after ~ 50 years it looks possible to translate between these two pictures, 
i.e., it seems possible to relate objects from these two perspectives. 

Therefore, a deeper investigation into the properties and features of QFT seems 
a reasonable thing to do. Studies in dualities (AdS/CFT, Langlands), emergent 
properties (effective theories), non-commutative structures, non-perturbative effects 
(branes), etc, seem like a rich venue to follow. 

The present work started as a means to deepen the understanding between symmetry 
breaking and the multiple solutions^ of a certain quantum field theory, in the spirit 
already developed in [35]. 

To better understand this line of argument, let us focus on 0-dimensional (i.e. 
ultralocal) models for the sake of simplicity. The idea is to compute the Schwinger- 
Dyson Equations (SDEs) of the theory in question and solve it analytically. Generically 
speaking, given that most theories are defined by a Potential which is a polynomial 
in the fields, the SDEs are expected to have more than one solution. Therefore, not 
all of these solutions will admit a series expansion in terms of the coupling constants 
(which are normally called "perturbative series" ) , implying that some solutions are 
non-perturbative, even though they admit a series expansion in some other parameter. 
In this sense, some of these solutions may keep the original symmetries (present in the 
Action), but some may not. To fix the idea, let us look at two particular examples. 



-"^ Henceforth, when referring to QFT and its solutions, we shall use interchangeably the following 
terms: solution space, configuration/phase space, vacuum manifold and moduli space. 
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1. Motivations and Introductory Remarl(s 



Airy Potential This theory is defined (in 0-dimensions) by >S'(0) = 0^/3+ J0, and 
its Feynman Path Integral (FPI) representation is given by, 

Z{J) = ^-^^ < oo ; 

e'^ d(t> 



c 

where C is the contour of integration that renders this FPI finite. On the other hand, 
its SDE is obtained in the following way: 

{4>^ -J)Z(J)^0 l^zlf^ (^2 _ J) 2;( J) = . 

Now we make an important observation: the SDE above is second order and, thus, 

must have two solutions. As it turns out, these solutions arc given by the Ai(J) and 
Bi(J) functions. At this point, it is not diflicult to realize that the FPI is nothing but 
the integral representation of the differential problem posed above, which implies that. 



Z[J] = = Ai[J] or Bi[J] 

JClo 



where the contours Ci,2 are such that the integral representation above is finite. 

The realization that should be accomplished at this stage is that 0-spacc had to 
be "complexified" so that we would be able to find both of these solutions — usually, 
when we think in terms of a FPI, it is presumed that the integration above should 
be done for e M. However, as we just showed above, in this very simple example 
we already sec that if we do not extend the FPI representation to G C, we will not 
obtain all of the possible solutions to the problem at hand. 

Having said that, it is important to realize that the partition function, in its most 
general form, will be given by a linear combination of all possible solutions, i.e., 

Z(J) = a Ai( J) + (3 Bi( J) ; 

where a, /3 e C are scalars. 

Further, the asymptotic series of Ai(J) and Bi(J) has a different form in different 
quadrants of the complex plane, a fact known as the Stokes phenomenon, implying that 
Z{J) will be a meromorphic function. Therefore, resummations of these asymptotic 
series, if it has to be done at all, must be tackled with the utmost care and attention, 
once mixing different branches of these series renders all results obtained after that 
completely useless. 

Also, because J G C, we can study the partition function in terms of its properties 
with respect to the Mobius group, which is the automorphism group of the Riemann 
sphere. In particular, we may focus in one of its more important subgroups, the 
Modular Group (SL(2,Z), where we identify the elements A and —A), and study all 
of the connections with modular forms, elliptic curves and many fractals — some 
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connections with String Theory can be drawn once it is realized that the generators of 
this group are given by "unit translation to the right" and "inversion in the unit circle 
followed by reflection about the line ^{z) = 0" , where this last generator is analogous 
to T-duality (which can be understood in terms of Mirror Symmetry). 

This is an interesting example because the cubic potential is usually said to be 
unbounded from below and, therefore, yields no meaningful theory. However, we 
clearly see above two solutions which are finite. This happens because we are explicitly 
looking for boundary conditions to solve the SDE, a fact that is completely analogous 
to searching for contours that render the Path Integral finite. 

Furthermore, this example can be readily extended from a scalar- valued field (f), 
to more interesting cases such as that of a matrix- valued field or of a Lie-algebra 
valued field ip. As expected, the only modification to the above equations is that the 
Airy function is extended analogously to either a matrix-valued or Lie-algebra- valued 
function — in fact, these extensions are defined via the FPI representation above [53]. 

Once these extensions are established, going from 0-dim to rf-dim is accomplished 
via the use of D-Modules: in this sense, the D-Modules are used to compute the 
particular "Thermodynamic Limit" arising from the 0-dim theory at hand. This will 
be important for future developments. 

It is important to emphasize that the observables of a theory should be real, once 
they are the measurable quantities. But this says nothing about the nature of the 
partition function nor of the fields themselves. 

Xcf)^ Potential This theory is defined by S{(f)) = ^i(f/2 + A0V4. The SDE and 
FPI for this model are: 

(/i + A 03 _ J) = (9^ + (7 9j - J) Z[ J] = ; 

J gi(M0V2+A</.V4-J</') ^0 

^[^] = ""'""r .2/2^, .4,,, = U[g, J] or V[g, J] or W[g, J] ; 

j e''^^'^ /^^^^ '^U<i) 

Cl,2,3 

where g = /i/A, and the three solutions for the SDE are given by the Parabolic 
Cylinder Functions, U[g, J], V[g, J] and W[g, J], which, in turn, can be obtained with 
three different contours, €1,2,3, in the FPI representation (see Appendix C). 

As mentioned earlier, each one of these parabolic cylinder functions represents one 
particular solution: symmetric, broken-symmetric and solitonic. 

Again, all of the considerations made above can be readily applied in this case: 
complexification of 0-space in order to obtain all possible solutions, the partition 
function is a linear combination of these three solutions. Stokes phenomenon is present 
in the asymptotic series for each one of these solutions, the partition function is 
meromorphic and, finally, its properties with respect to the Modular Group can be 
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1. Motivations and Introductory Remarl<s 



studied. 

1.1. The road ahead 

We have just shown how extremely important boundary conditions (resp. contours) 
are for properly defining a QFT and all of its solutions/sectors (resp. moduli space); 
and this is the cornerstone of our work: from this point on, we can attack several 
different problems in novel ways. 

Here are the problems that we tackle in this thesis: 

Mollifying QFT We probe these multiple solutions using the help of a mollifier, which 
is essentially a band-pass filter that eliminates the meaningless highly-oscillatory 
parts of the Partition Function, but does retain its meaningful structure. This 
shows clearly the links with Lee- Yang Zeros (complex zeros of the Partition 
Function, seen as a polynomial in the coupling constants) and Stokes phenomena 
(the asymptotic expansion of each solution has a parameter-dependent limit, 
which means that there is a discontinuous change [in the Partition Function] 
when crossing Stokes lines separating distinct solutions/sectors of the theory). 
For more on this, please see [40]. 

Topological Transitions and Geometric Langlands Duality We show that each dif- 
ferent solution/sector of the theory has a distinct topology which can be used as 
a Superselection Rule (topological charge); moreover, the Euler Characteristic 
(second Chern class) of each phase gives its quantization condition. Also, these 
sectors can be labeled by the discriminant of a polynomial constructed from 
the potential, where the structure of this discriminant is related to singularities 
in the Partition Function, which, in turn, must be treated as a meromorphic 
function (which connects with the picture above about Lee- Yang zeros and 
Stokes Phenomena). However, the different solutions of the same theory are 
related by "dualities" in the values of the parameters of the theory: distinct 
boundary conditions establish different allowed ranges for these parameters. In 
this sense it is possible to construct several "dualities" that translate one sector 
into another, and this is how we make touch with the subject of Geometric 
Langlands Duality. Furthermore, using Z)-Modules, we can take the Thermo- 
dynamic Limit of each particular 0-dim theory into its full d-dim version, thus 
dimensionally constructing each different sector of the theory: each branch is 
"grown" out of the "natural differential operator" defined by the specific poten- 
tial and its parameters (determined by the appropriate boundary conditions) . 
Therefore, these dualities are not only with respect to different sectors of the 
theory, but also with respect to each natural differential operator associated 
to each phase, thus establishing the possible allowed values for the geometric 
quantities (connection) associated to each point in the moduli space of the whole 
theory. For more on this, please refer to [41, 55]. 

Three-dim Gravity and Phase Transitions This is where we make use of the exten- 
sion to Lie-algebra- valued fields mentioned above when solving the Airy potential. 



1.1. The road ahead 
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Recently, Witten has shown that 3-dim gravity with negative cosmological con- 
stant has a dual 2-dim CFT, [39]. In order to do so and obtain a holomorphically 
factorized Partition Function, he had to move from real saddle points (which 
do not account for the whole theory) to complex ones. And, in doing so, he 
found that the Hawking-Page phase transition is analogous to the condensation 
of Lee- Yang zeros. The question we are tackling is that of extending the Airy 
function as previously discussed (to a Lie-algebra valued Airy function) and 
"grow" this 0-dim solution into a 3-dim one using its D-Module. Furthermore, 
there are some questions that arise from the connection with Braid and Knot 
theory, such as whether having two distinct Partition Functions (for Ai and Bi) 
imply different polynomial invariants, or whether this implies that the structure 
constants of the Lie group in question can be dynamically determined (by the 
allowed contours), etc. 



1. Motivations and Introductory Remarl(s 



Part I. 

Mollifying Quantum Field Theory 
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Chapter 



2 



Motivation and Introduction 



2.1. Sign Problem 



One of the fundamental difficulties of Monte Carlo (MC) approaches is known as 
the "sign problem" . It is encountered when the functional integrals to be evaluated 
do not have a positive definite measure. It is not related to any approximations 
or fundamental errors in the MC scheme, but it describes the situation where the 
statistical error can become very large. In general, any expectation value can be 
written as. 



where and O are real-valued functions of the field variables and, in general, the 
measure need not be positive: in Lorentzian QFT this measure is complex- valued 
and given by e^'^'^^^D^. However, if changes sign, it cannot be considered a 
probability density (z.e., a measure). 

The standard trick to avoid this problem is to modify the measure in the following 
way: = |/i[0]|/ J |/^[0]|- Then, one absorbs the sign of in the quantity to be 
measured: 



In some cases this may work. However, the random walk guided by is very likely 

to predominantly sample unimportant regions in phase space. 

At this time no completely satisfactory solution to the sigh problem exists, al- 
though there are some promising attempts including positive projection [ , >], fractal 
decomposition scheme [1], and Berry's phase ["] and its Stiefel manifold [fi]. We 
believe that the techniques discussed in the present work lay the grounds for another 
possible solution path. The general success of the methods studied here will depend 




(2.1) 




(2.2) 



9 



10 



2. Motivation and Introduction 



on the development of more computationally efficient algorithms than the ones used 
to illustrate the simple examples treated below. 

2.2. Smoothing out the Measure 

The idea behind the mollification technique is to use a convolution in order to smooth 
out and filter the measure: the highly oscillatory measure given IS 
convoluted with some suitable function (called a mollifier; see Appendix A for more 
details) and as a result only an effective contribution is left (rapid oscillations of the 
measure cancel out when integrated over the slowly varying mollifier). 

A convolution is an integral that expresses the amount of overlap of one function 
g as it is shifted over another function /. It therefore "blends" one function with 
another. In mathematics, moUifiers are smooth functions with special properties, used 
in distribution theory (generalized functions) to create a sequence of smooth functions 
approximating non-smooth functions, via a convolution. 

In an application of this idea to QFT we will replace the rapidly oscillating measure 
in Minkowski space (/ = exp{i S[(j)]}) by smooth functions (/^ = r]^* I) so that we can 
run numerical simulations that are otherwise prohibitive (see A). At this stage, the 
only requirement that is made is that: J r]^{x) dx = 1, where the integral is performed 
over the domain of rj^. The parameter e controls the approximation, the smoothness, 
being made, and in the e — * limit, the original expressions are recovered. (The 
meaning of this e parameter is that of dilating or contracting the mollifier itself, thus 
controlling the range over which the filtering is done. Note that the calculations should 
be done using a non- vanishing value for e in order to avoid the analytical expressions 
which are hard to handle numerically.) 
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Mollifying Quantum Field Theory 



The technique described above (and detailed in A) can be apphed to [Lattice] QFT in 
the fohowing way: 



This form is useful to set up numerical computations. Performing the mollification 
before performing the path integral can be very advantageous. 

Mollifying the integrand (i.e., the complex exponential of the action) smooths the 
highly oscillatory integral. The integrand changes its form, from the canonical / to 
the mollified Is, in the following fashion. 



i.e., the convolution with r]^ changes variables: ip \^ (p. Taking [functional] derivatives 
of Z (which yield Green's functions) is just the same as taking them with respect to 
because the derivative operator commutes with the mollification. 
The [Feynman] Path Integral is constructed just as before: 




(3.1) 



J[0; J] =e^S[0;J] 



J] = iv, * I)[if; J]= ^r],[ip- 0] e^^l-^^^l 2)0 , 



Z[J] = J/[0; J] 2)0; 
ZJJ] = J/,[v9; J]2)(^ 



Note that. 




hi-fi; J] 



.-. Z,[J] = Z[J] ■ 



(3.2) 
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3. Mollifying Quantum Field Theory 



where standard properties of convolutions (see [ ]) have been used in the last step. 

Even though the above result is analytic, when one goes to the simulations, a small 
dependence on e shows up. (See Section 4 for a more detailed discussion on this 
matter.) 



3.1. Importance Sampling 



The next step consists of choosing an appropriate sampling function (for the Monte 
Carlo simulation): 



where TiW = W[lp] "Dip, i.e., just a change of variables in order to make the computa- 
tions more robust: it is possible to use the integrand profile in order to speed up the 
calculation. 

An appropriate choice for the importance sampling function, W, is usually given by 

A little digression is in order: If all we wanted to do was to simulate a certain 
[Euclidean] QFT on a lattice, the above reasoning would be just as valid, modulo the 
mollification, i.e., we would conclude that the proper importance sampling function 
was W[(p] = I / [0] I • The full-fledged formula is useless for doing computations since 
it implies full knowledge of the theory being calculated. A better choice is a simple 
approximation: we choose a saddle-point approximation in order to make the expression 
more manageable and still keep most of the characteristics of the integrand. 

Thus, the importance sampling function (Ty[0]), which is used just to better guide 
the simulation, would be nothing but the absolute value of the sum over all saddle- 
points of the theory in question. Assuming we have only one such saddle-point, the 
well known answer is given by: 



W\ 



3* SIM 



det{-92 + m2 + V"[0o]} 

-itr{log(-92+m2+V"[</.o])} 



e 2 



|tr{log(-a2+m2+V"[9io])} 



det{-92 + m2 + V"[0o]} 



-1/2 



(3.3) 



Using this to find an importance sampling function for a simple scalar QFT on the 
lattice requires solving for the determinant above. In Lattice QCD, the problem of 
the fermion determinant is a very time-consuming operation that must be done for 
every step of the MC calculation. 

But this is not the only computational bottleneck. The more severe one comes 
from the so called sign problem: The exponent in is not bounded from below, 
therefore one cannot guarantee the ergodicity of the Markov Chain underlying the 



3.2. Mollifying the Importance Sampling Function 
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Monte Carlo draws, i.e., a much bigger number of draws would be needed in order to 
yield any meaningful answer. This is the reason why calculations are done in Euclidean 
space rather than Lorentzian/Minkowski space. In Euclidean space (Wick-rotating 
the integrand), the exponent in question becomes = e""^'"^!; it is bounded from 
below and the Monte Carlo method works fine (the Markov Chain behind it becomes 
ergodic). 

It is to tackle this problem that the method of moUifiers comes in: mollifying 
will yield a smooth function, whose properties help the convergence of the MC 
computation. Furthermore, if this can really be accomplished and calculations with 
imaginary exponents become a reality, the next logical step is to analyze the different 
phases of the given QFT. As shown in [0], the different phases of a theory can be 
picked out with a mere choice of boundary conditions, which is the same as properly 
defining the measure of the [Feynman] Path Integral. The upshot is that if we want to 
obtain answers other than the ones that can be reached via perturbation theory, either 
the measure of the path integral or the boundary conditions of the Schwinger-Dyson 
equations have to be modified [')]. 

Once there is nothing preventing the measure of the [Feynman] Path Integral, from 
being (in the most general case) complex-valued — in order to account for the phase 
structure of the theory — , the fact that its integrand is also complex- valued has to be 
taken more seriously, afterall this is a non-perturbative result: a perturbative series 
only works if we know, a priori, in which phase we are working, so it can be tailored 
to that particular sector of the theory. 

Thus, the need to address the sign problem becomes even clearer: it is not simply a 
problem of doing MC in Lorentzian/Minkowski space, it is mostly a problem of being 
able to compute all possible solutions — all the different phases — of a given QFT. 

3.2. Mollifying the Importance Sampling Function 

As outlined before, the idea is to use the mollification technique to handle the highly 
oscillatory terms present in the importance sampling function. 
There are three possible choices for the sampling function W: 

1. Taylor-expand and mollify the complex integrand: no knowledge about the 
saddle-points is necessary, but it is very computer intensive, which is the reason 
we will not focus on it; 

2. As done in (3.5), where the integrand is saddle-point expanded and mollified, but 
the explicit form of the mollifier is used in order to perform the remaining integral: 
not so heavy on the computer, but some knowledge about the saddle-points is 
needed; 

3. As done in (3.6), where the integrand is saddle-point expanded and mollified, 
generalizing the standard saddle-point approximation: less intensive of all 3, but 
even more knowledge about the saddle-points is needed. 
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3. Mollifying Quantum Field Theory 



The difference between the second and third methods above is that the integral in 
(3.4) is exphcitly carried out in the former, but not in the latter. As explained below, 
the third method (listed above) consists of the simple generalization of the well known 
saddle-point approximation; however, we need to note that more knowledge about the 
particular saddle-points is needed in order to perform the contour integration involved 
in this scheme (and this is the difference with respect to the second method). 

The calculation below shows the second method outlined above, known as 2 "''-order 
gradient approximation: 



Now, a saddle-point expansion is performed on the action. 



A moUifier needs to be chosen, and for the purposes of this calculation, a Gaussian 
one is our choice: 



-oo 

exp \i S[ipo] + ('^-'^o)^"[^o](0-^o) 



V27re2 



2 l-ieS"[^o\t 



det{l-ze5"[v9o]e} 



(3.4) 
(3.5) 



The important point to note about the derivation above is that the explicit form 
of the mollifier, 77^, had to be used and, because of that and the nature of the MC 
simulation, we do not necessarily need a good knowledge about (pQ. This means that 
we can compute the above importance sampling function, (3.5), and use a trial (/Pq''"^^: 
the [MC] simulation will do the job of moving towards the exact ipQ and pick out the 
different phases of the theory. As shown elsewhere, [11, 12, 13, 8], this works quite 
fine. 

However, for more relevant cases (e.g. 4-dimensional QFTs) the time required for 
simulations becomes unrealistic. The bottom-line is that a simple-minded implemen- 
tation of the method will have the code/computer doing the work of mapping the 
(0, v9)-space, while one could use that information beforehand in order to speed things. 

The way out of this is to generalize the following well known result: 3{s) = 
J^g{z)e'f^'Uz = g{zo) e' ^'^''^'^ e'" / f"{zo)\ -, seR,zeC; i.e.. 



iy,[0] 



'det{5"[^o]}| 



(3.6) 



3.2. Mollifying the Importance Sampling Function 
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At this point, the similarity between (3.3) and (3.6) is clear, in fact, they are the 
same formula, except that, in the former case there is no mollification while in the 
latter there is mollification. 
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Chapter T 

Tuning the Mollification 



At this point we are almost done with the analysis of the mollification technique, the 
only aspect remaining being the optimal choice of e, the parameter that regulates the 
approximation. 

From the discussion in appendix A, it is clear that in the e — > limit we recover 
the original theory. However, numerically, this is not exactly what happens, because 
of numerical and statistical fluctuations. Therefore, the choice of an optimal value for 
e is crucial. 

Using an information theoretic viewpoint, we can choose e to optimally compress 
^^£[0] around each stationary phase point. Therefore, we can define an information 
entropy based on in the following way: 



J (4.1) 



The optimal e is the one that minimizes S[e]. This value corresponds to maximal 
compression of the information in the importance sampling function. 

In order to illustrate the above, let us consider a very simple [0-dimensional] example, 
given by S[(j)] = 171"^ (p"^ /2. In this case, the saddle-point is 00 = and the importance 
sampling function and the information theoretic entropy are given by: 



S[e] 



m e v2 TT 
log(m e V2 TT ) — 
1 — m log(m) 



Plotted below are the graphs of the importance sampling function (also showing 
how it varies with e) and of the information theoretic entropy (both using m = 1): 
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4. Tuning the Mollification 



IrrpDrtinofl aampirtg rundlnn: W^^il 




a a.4 a.E d.b 



Figure 1: Importance Sampling function (leftmost plot, m = 1, Wc[0] = e /^"^ /e tt ) and 
information theoretic entropy (rightmost plot, m = 1, S[e] = log(e \/2 tt ) — e^) for the action 

= 02/2. 

We can clearly see the importance function peaking around (j) = 0, which is the 
expected behavior {i.e., it happens around the saddle-point 0o = 0), while the 
information theoretic entropy shows that the best parameter for the mollification is 
e — > 0: this is not desirable, once lime_+oS(e) ~^ —oo, but it simply shows that we 
should use a small non-vanishing value for e, as mentioned before in Section 2. 



Chapter 

Simulation Details 



In order to perform the simulation as outlined in 5.4, some technical details had to 
be taken into more serious account. Among them, the important ones are: random- 
number generation and non-local algorithms. Those will be the bottle-necks, as noted 
in [11]. 

5.1. Random Number Generators 

Even though some consider the random-number question a solved one, this is not 
always the case [ ]. For this reason, the choice of which random-number generator 
(RNG) to use is still a critical one, not only for speed reasons but because of systematic 
errors as well. In order to address this issue the choice made was for the Mersenne 
Twister RNG [15] (period of 2^^937 _ 

5.2. Non-local Algorithms 

As for non-local algorithms, the choice made was for a Genetic Algorithm (GA) [16]. 
The "population", in the Lattice QFT sense, consists of the lattice {i.e. spacetime) 
points. For a given "recombination rate" , an initial population gets randomly divided 
into pairs of [lattice] points and "crossing-over" (or "genetic" ) operators are applied to 
these pairs in order to generate different sets of pairs. In this context, the crossing-over 
operators exchange the coordinates of given points, thus there are 4 of those operators, 
namely: T^, T^, Ty,Tz. That is, if the crossing-over operator is applied to the 
points (p = {(f)t,4>x,(py,(pz) and (j)' = (f)'^, (j)y, (p'^) the outcome will be: Tx[4>,(f)'] = 
{{(pt, (p'x, (py,(pz), (0t, (px, 4>y^ (f^'z)}- Analogous definitions are valid for the other crossing- 
over operators. The key features for this choice were: which Tj to use is random; 
whether or not the change i— > is made depends on a given probability and the 
recombination rate can be arbitrarily chosen (although keeping it below 1% showed 
to be a good tune). Moreover, the genetic operators are unitary, i.e., Tj = Tf^. 
This implies that = 1, which guarantees the so-called detailed balance of the 
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5. Simulation Details 



MC simulation (this is the ergodicity of the algorithm). At this point, the only step 
remaining is the explanation of the probability profile used in algorithm 2 below. 



1 


Choose the recombination rate such that: ^ r ^ 1 {usually between 0.5% and 




1.0%} 


2 


Draw a random number /5 G [0, 1] and compare with r 


3 


if /3 ^ r then 


4 


Metropolis Monte Carlo loop {using the random walk technique} 


5 


else 


6 


Genetic loop {see algorithm 2 below} 


7 


end if 


Algorithm 1: Genetic Algorithm 


1 


Choose a random pair: {0i, ^2} = {0(^1, yi, zi), 0(^2, X2, ^2, -22)} 


2 


Generate the pair {0'i,02} = {(l){t[, x[,y[, z[), 0(t'2, I/2) 4)} = ^i(0i,02) 


3 


Draw a random number c G [0, 1] 


4 


if c>P(0;,0'2)/P(0i,02) then 


5 


do nothing 


6 


else 


7 


perform the exchange (4>i,(f>2) ^ (0'i)02) 


8 


end if 



Algorithm 2: Genetic Loop 



In the notation used below, P is such probability profile and P{4>'i, (p'j) means that 
this profile is calculated using the [GA generated] points 0^ and 0^. It proved useful 
to chose the mollified importance sampling function as this profile P. Also, note that 
0j and just like 0'j and 0^, are just particular values of the field at the sites i 
and j. A comparison between these two profiles is then performed. 



5.3. Euclidean Lattice Simulations 

In order to be able to understand better what will happen in the Lorentzian/Minkowski 
case, a quick and dirty reminder of the Euclidean one is presented below. 

In what follows, $ is the set of all possible field configurations and H^H is its 
cardinality (i.e., the number of elements in the set). For each element of $ (i.e., for 
each field configuration) one has to average the observable over all the lattice points 
(which are computed via Metropolis). Each element in the Markov Chain (generated 
by the Metropolis MC algorithm) is denoted by 01*'. To obtain the final result, an 
average over all field configurations is made. (To make clear that the final goal is the 
actual lattice computation, the measure is denoted as [c/0].) 



5.4. Mollified Minkowski Lattice Simulations 
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JO[0] exp{-5[0]}[#] 
Jexp{-5[0]} 

1 
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Volume 




Volume 
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5.4. Mollified Minkowski Lattice Simulations 

In the Minkowski version of the above, the Wick rotation is not performed, therefore 
we are left with the original form of the exponent, 

JO[0] exp{^^[0]}[#] 
Jexp{i^[0]}[#] 

The above functional will be moUified in order to yield more tractable expressions. 
The fact that there will be two functional integrations (rather than just one, like 
above) should not bring many problems. The real question here stems from the fact 
that the moUifier, 77, mixes the two of them. 



\[\ne[ct> - OM exp{i5[(^]} m][dcl>] 
J{J?7e[0-<^] exp{i5[(^]} [o?(^]}[#] 

ll^ll ll^ll h I S£=i ^^[^^^ - exp{z5bW]}/W^f M']] J ' 

(5.1) 

where it is understood that the configurations in $ are chosen with respect to the 
importance sampling function, and the configurations in F are chosen with 

respect to a uniform distribution. This happens because the mollification process 
(namely the convolution) mixes the variables from $ and F together. Therefore, the 
only way to implement this "mixing" is by having a uniform distribution for F and 
implementing the "interaction" with the variables in $ — via the mollification — 
explicitly, using 77. Note that ||$|| = \\F\\. 



{0) = 
•■• (O) = 
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Chapter 



6 



Spontaneous Symmetry Breaking 



We have been discussing lattice QFT in Lorentzian/Minkowski spacetime because one 
of the main objectives of this paper is to examine a numerical approach to directly 
calculate Green's functions in different phases in QFT. This necessitates being able to 
evaluate path integrals of complex exponentials. The arguments for this have been 
given elsewhere [9] and will be illustrated in a particular example in the following 
section, but it is fairly easy to understand why this is the case. The traditional 
[Feynman] Path Integral formulation of a QFT involves integration of the exponential 
of the action over every field variable at every space time point. Traditionally, these 
integrations (assume the action is written in terms of self adjoint fields) range from 
negative to positive infinity along the real axis. In the limit of small couplings, this 
form of the [Feynman] Path Integral generates perturbation theory and the results 
appear (for finite number of expansion terms) to be regular at vanishing coupling. 
Thus, all expansions that are not regular at vanishing coupling, such as the traditional 
symmetry breaking expansion of quartic scalar field couplings, are excluded. In order 
to avoid this restriction and produce all possible solutions of the QFT it is necessary 
to extend the [Feynman] Path Integral integrations to complex values of the fields in 
a way consistent with the field equations and reality properties of the theory. 

Indeed, as shown in ['!] and section 7, the different phases of the theory emerge 
through the varied boundary conditions consistent with the equations of motion. 
Equivalently, rather than varying the boundary conditions of equations of motion, 
the measure of the Path Integral can be changed. In general the number of choice of 
paths of integration for the path integral correspond to the number of independent 
solutions of the Schwinger-Dyson differential equations. 

If we define a QFT via its [Feynman] Path Integral, all we need to know is the 
action, S[(f), J] (for a field whose source is J), of a given model, for then we can 
write the generating functional as: 



where !N is a normalization constant such that Z[J = 0] = 1. The crucial question that 
remains unanswered in this approach is: "How does one properly define the measure 
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6. Spontaneous Symmetry Breaking 



V(f)?" 

The best answer so far (for 4-dimensional systems) says that this can only be done 
for free QFTs, via the use of cyhndrical functions [17, 18]. (Note that this only happens 
in the continuum. In its Lattice formulation, QFT is free from such peculiarities 
because of the lattice regularization.) 

Analogously, defining a QFT via its Schwinger-Dyson equation, 

_ j(^) z[J] = ; 

is equivalent to substituting (J) —i jj in the action, requiring the equations of motion 
to be the solutions that extremize it. Note that the Schwinger-Dyson equations are 
a system of infinitely many [partial] differential equations, one per each point of 
spacetime. 

Therefore, when thinking in terms of differential equations, the boundary conditions 
are responsible for the phase structure of the theory. Intuitively, the picture that 
comes to mind is that of a portion of space divided into as many subsets as there are 
solutions to our equations of motion, such that in each of those regions, the equations 
of motion satisfy appropriate boundary conditions. 

On the other hand, when thinking about the [Feynman] Path Integral, usually it 
does not seem bothersome that the measure is not properly well-defined. In fact, 
determining the boundary conditions for the Schwinger-Dyson equations is analogous 
to determining the measure for the Path Integral. Thus, just like the boundary 
conditions, the measure is responsible for the phase structure of the theory (in the 
integral representation of the problem). 

For completeness sakes, this is how a QFT in different spacetime dimensions, d, 
behaves: 



d = For 0-dimensional QFTs — i.e., QFT on a point (the [Feynman] Path Integral 
degenerates into a simple integral) — there is no such thing as phase transition 
since there is no such thing as dynamics. However, the phase structure of the 
theory survives, given by the different boundary conditions (resp. measure) 
needed in order to determine all the solutions to the equations of motion. 

d = 1 For 1-dimensional QFTs — i.e.. Quantum Mechanics — , again, there is no 
such thing as phase transitions, for if the theory has 2 different vacua we could 
take a linear combination of them to be the "real" vacuum state, given that they 
would be related by tunneling. (For the analogous case in Condensed Matter 
physics, please refer to [if].) 

d ^ 2 In this case, phase structure and transition exist; this is the complete scenario. 
(Note that there is no operator that relates 2 inequivalent 6'- vacua of the theory, 
therefore they belong to different algebras.) 



The real lesson to be learned from all of this is that Spontaneous Symmetry Breaking 
is a phenomena generated by boundary conditions; whether you use them to define 
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the measure of the Path Integral or to define the Schwinger-Dyson equation is just 
a matter of personal preference. It is not said anywhere that the limits of the Path 
Integral have to be real; what we need to have are real observables. In fact, as shown 
below on section 7, it turns out that in order to have symmetry breaking we need 
a measure that is not necessarily real: this will enable the computation of all the 
solutions of a given QFT [9]. 
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Chapter I 

Results 



In what follows, the results obtained thus far are presented. In short, they are in 0, 1 
(time) and 4 spacetime dimensions. 



7.1. Lower Dimensional 

Let us start by addressing the 0- and 1-dimensional results. Plainly and simply put, 
this means that we are solving a simple integral in spacetime dimensions since the 
[Feynman] Path Integral degenerates into a standard integral and, in the 1-dimensional 
(time) case, one will be doing Quantum Mechanics. 

For the 1-dimensional results, refer to [8]. There, the Quantum Chemistry of the 
problem is fully treated and addressed. Note, however, that the importance sampling 
function chosen in [ ] is different than the one used in this work: the quantum 
chemistry was done using (3.5) while we use (3.6). 

Below, the 0-dimensional results are summarized and, in order to illustrate the 
features of the moUifier technique, different properties of those models are made 
explicit. 



7.1.1. Airy Function 

The action for this model is given by: S{x) = + J x. This is an interesting model 
because we can explicitly calculate the partition function and compare it with the 
results coming from the Mollified Monte Carlo procedure. Moreover, this model has 
2 stationary phase points and the results displayed are from the one in the complex 
plane which is not accessible with normal Monte Carlo. The quantity of interest is: 
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AkirFDUrtldl 




Figure 2: Plot of the Airy potential, V{x) = highlighting the fact that it is not bounded from 

below. 



Z[J] 



/ expo \-iJx>ax 



Ai(J) 



|_^exp{zi^}cix 



, Ai(0) 



The [first two] graphs below show the highly oscillatory behavior of the integrand 
in the partition function: the left one is its the real part, while the right one is its the 
imaginary part. 




Figure 3: Real and Imaginary parts of q^^'^/^+^J^ showing the highly oscillatory nature of the 

problem. 

Below we see the graphs of a particular 2-dimensional slice of the above pair, where 
J = —1: the real (left) and imaginary (right) part of the mollified Airy-integrand (in 
red, e = 0.1) is in contrast to the non-mollified integrand (green). 



7.1. Lower Dimensional 
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Figure 4: Superimposed plots of mollified (red) and non-mollified (green) real and imaginary parts 
of the integrand e"^/'^"", showing the smoothness achieved on the J — —1 two-dimensional slice of 

the previous graphs. 

Below we see the graphs of the real (left) and imaginary (right) part of the mollified 
Airy- integrand for different values of the field x and of the mollification parameter e: 
analogous to Figure 3 but for the mollified version of the integrand (which corresponds 
to the red lined plots of Figure 4). 




Figure 5: Plots of the real and imaginary parts of the mollified integrand for various values of x 

and e. 

The 2 stationary phase points are given by xq G {±^/— J }, where xq is the solution 
of S'{xo) = 0, and the [properly mollified] importance sampling functions are: 

W^/^°N=exp{-^;^}+exp{-l^}; J^O; (7.1) 

wn^] - -p{-l (^^)] +exp{4 (^^^)] ; ^ < ; 

(7.2) 

It is customary to see the source term above called t and identified with time but this 
is not the approach taken here. 

Below we see the graphs of with respect to the e parameter and x, for fixed 
values of the source: on the left J = 4, and on the right J = —16, 
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Figure 6: Plots of (7.1), for J = 4, and (7.2), for J = -16. 
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Figure 7: Plots of 2-dimensional slices of (7.1) and (7.2) for various values of e. 



It is not difficult to see that there is an optimal value for the parameter e, as shown 
below in the graph for the information theoretic entropy (see Appendix B for further 
details): we have the information theoretic entropy for J ^ on the left, and J < 
on the right, 



Iniyxmalkin TtiDcnilt Enlrcpf: 3 



Inrcrmalkn TMDDffflt ErUrEf^ E i^l 











































































— ' 


\ 












Figure 8: Information theoretical entropy for finding the optimal value of the mollification 
parameter: on the left we have J ^ 0, and on the right we have J < 0. 



7.1. Lower Dimensional 
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Furthermore, it is easy to perform the analytical calculations outlined above (using 
the Airy functions on the Generating Functional, rather than its mollification) to 
get the pure evaluation of this model. The results show the mollified answers (given 
a proper choice of e) are sensationally accurate and show, already in this simple 
example, that a region not allowable in ordinary Monte Carlo approaches can be 
computed. In fact, due to the simplicity of this example, modern computers can 
calculate it straightforwardly, without the need of special tricks to handle the highly 
oscillatory integrand. The point of this example is twofold: show the precision and 
the accuracy that the moUifier method can achieve and also to illustrate how this 
smoothing procedure works. 

7.1.2. 0-dimensional 0^ Theory 

Also known as ultra-local 0^, its action is given by: S[(f)] = /i 0^/2 + (? 0^/4. As before, 
the saddle-point 0o is such that S"[0o] = 0. The graphs for the Action above (positive 
mass on the left, negative mass on the right) are given by: 




Figure 9: 0^ potential for positive and negative values of fj,. 

Below, we show the highly oscillatory behavior of the integrand of the Partition 
Function given by the above action, i.e., / = e*^'^'^ /2+5</' /4)^ ^^j. positive and negative 
values of the fi parameter: 
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Figure 10: Plots of the Real and Imaginary parts of the integrand for various values of J and 

^ — 1, g — I and fi = —1, g = 1. 

Below, we see the real (left) and imaginary (right) parts of the mollified 0^-integrand 
(red) in comparison with their non-mollified counterparts (green): 
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Figure 11: Superimposed mollified (red) and non- mollified (green) real and imaginary parts of the 

integrand for fi ~ I, g = I, J — —1, e — 0.1. 

Below, we see the real (left) and imaginary (right) parts of the mollified 0^-integrand 
for different values of the field and of the mollification parameter e: 
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Figure 12: Plots of the real and imaginary parts of the mollified integrand for various values of the 

field and the parameter e. 



7.1. Lower Dimensional 
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As explained in [9, 12], the three different solutions/phases of this model can be 
selected via an appropriate choice of boundary conditions: = M, F"*" = (— oo, 0) U 
{0,ioo) and F~ = (— oo,0) U (0,-200). The leftmost figure below shows those 
boundaries in the Argand Plane: F*^ = M is the green line, F+ = (—00, 0) U (0, i 00) 
is the red line while F~ = (—00, 0) U (0, —i 00) is the blue line. The rightmost figure 
below shows the regions (shaded) of the complex 0-plane, (0 = pe*^), defined by 
cos(4^) ^ 0. Any contour, starting and ending at infinity, within one of these four 
domains corresponds to a particular solution of the 0-dimensional (j)"^. 




Figure 13: Left: Three different contours thiat render tire Patlr Integral finite and depict the three 
distinct solutions obtained in this model. Right: Any contour that starts and finishes inside the 
same-color shaded regions (following the general directions shown on the left) corresponds to a 

particular finite solution. 



Using the three different contours above we are able to find the three different 
solutions to the equation of motion of this 0-dimensional QFT. The Schwinger- Dyson 
equation is given by, 

-^(^^- + g—j Z[J] = JZ[J], 

where we can easily see (again) that there has to be three solutions. These are 
called Parabolic Cylinder Functions and are denoted U{fi/g, J), V{fi/g, J), W{fj./g, J). 
Those behave in the following manner (see Appendix C): 

Regular at g — ^ Consistent with perturbation theory. 
Singular (2± ^/g^) at P — Symmetry Breaking, 
Singular (c± exp{fi/4g}) at g — ^ Instanton. 
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Thus, 

oo 



mj] = 



-oo ro.± 



It should be noted, though, that depending on the chosen contour (r°,r+,r^), 
the above integral will resemble a Presnel one and, as such, its computation is quite 
simplified. 

For the graphs below, we used the following values for the saddle-points, 
where <S"[0o] — 0. Therefore, the difi^erent importance sampling functions are: 



Symmetric Phase Given by the r° contour, this is the one which is regular in the 
limit g ^ and accessible via perturbation theory, 

<[0]=exp|-i(^)'|; 00 = 0; 



Solitonic Phase This phase is given by the linear combination of the contours F"*" 
and F~ such that ^ 0, i.e., this represents the soliton solution. 



Broken-Symmetric Phase This one is given by the linear combination of the contours 
F+ and F~ such that < 0, i.e., this represents the solution usually referred to 
as broken-symmetric, 

WrVl=exp{-i(*±^)^}+exp{-i(i^)'}; ^- = ±^=7^; < 0) 



Below, we have the graphs of these importance sampling functions: 



7.2. Higher Dimensional 
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Figure 14: Plots of the mollified importance sampling functions for various values of the field and e: 
symmetric ((/)o — 0), solitonic (fx — I, g — 1) and broken-symmetric (/j, = —1, g — 1). 

The graphs below show the entropy for the above importance samphng functions, 
as discussed in Appendix B: 




Figure 15: Information theoretical entropy for finding the optimal e in the symmetric, solitonic and 

broken-symmetric cases. 

7.2. Higher Dimensional 

As for the 4-dimensional results [for a Free Scalar QFT] the news is not encouraging. As 
mentioned in section 3.2, there are severe efficiency constraints in this implementation 
of the mollification; in section 3.2 three possible choices of solutions to this technical 
issue were given. 

It turns out that, from those three solutions, one is non-efficient while the other 
two produce analogous results: only the imaginary parts of physically interesting 
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quantities (generating functional and Green's functions) can be tamed. The real parts 
continue to show an oscillatory pattern. 

Although this is at least disappointing, the solutions obtained are good enough to 
control the imaginary part's oscillations. 

It is our hope that changes in algorithms will allow some version of mollification to 
work in real world problems. 



Part II. 

Solution Space and Moduli Space 

Topology 



37 



Chapter v>l 

Introduction 



Classical gauge theories are well described through differential geometry, where a 
gauge field is represented by a connection on a principal fibre bundle P for which the 
structure group is the symmetry group of the theory; see, for example, [22, 23, 24, 25, 
26, 27, 28, 29, 30] and references therein. 

The phenomenon of symmetry breaking also has its own geometric formulation 
[25, 26, 27, 28, 29, 30, 31]: the reduction of the principal bundle P. To this classical 
picture, the question that arises is that of how does a quantum field theory select 
its vacua given by the different possible reductions of the principal bundle P. The 
canonical answer to this question is that radiative corrections (beyond the tree-level 
approximation) cause the QFT to "jump" from the symmetric phase to the broken- 
symmetric one(s) [ ]. 

However, [ ; ;] already expresses concerns about some issues (e.g., the last paragraph 
on page 1894 and its continuation on page 1895) that are more explicitly treated 
in [35]. For more examples about these issues, see [o()] (pages 27 and 28 discuss 
the asymptotic nature of the series expansion, their maximum accuracy and issues 
regarding Borel summation) and [12]. 

Loosely speaking, what [35] does, is to compute the Schwinger-Dyson equations (for 
a given QFT) and note that these differential equations have as many solutions as 
its order indicates. Moreover, the boundary conditions that determine each of these 
solutions yield different possible values for the parameters of the theory (mass, coupling 
constants, etc). Thus, each one of these solutions has its own series expansion (which 
may or may not be equivalent to the perturbative series) and particular behavior. 

The analysis of the boundary conditions of the Schwinger-Dyson equation being 
responsible for the different solutions of the theory done in [35] is more in the spirit of 
the self-adjoint extension of the associated Hamiltonian, as done in [13, 44]. Although 
it may not seem so at first, [ ; ] is equivalent to a latticized approach to QFT and, as 
such, requires 2 types of infinite limits in order to give rise to its continuum version: 
the thermodynamic (number of particles in the box) and the volume (size of the 
box) limit. These limits do not commute and fiddling with them (as done in [35]) is 
analogous to resumming the perturbative series (as done in [ ' •]). 

The present work has the goal of developing a geometrical method that brings to the 
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8. Introduction 



foreground these issues of vacuum structure in QFT, its multiple solutions, its moduli 
space, clearly showing that these different solutions are topologically inequivalent. 
The fact that different configuration spaces for distinct solutions of the equation of 
motion of a given Lagrangian have different topologies shows that expansions must be 
performed separately for each solution of a theory, i.e., each phase has to be regarded 
and treated as a separate theory. 

In this sense, this work streamlines how the parameters of a theory (mass, coupling 
constants, etc) determine the topology of the vacuum manifold (moduli space) which, 
combined with the picture presented in [35] , gives a prescription of how to examine 
the solution space of a field theory: The boundary conditions of the Schwinger- 
Dyson equations determine the parameters of the theory which, in turn, determine its 
topology. 



Chapter 



9 



Classical Field Theory and the Jacobi 



Metric 



A typical action for a scalar field has the form, 



= J K{n,7T)-Vr{<f>)(rx; 



where K is the kinetic quadratic form (bilinear and, in case tt e C, hermitian; defining 
the inner product g) and Vt is the potential, where the index r collectively denotes 
coupling constants, mass terms, etc. 

We want to reparameterize our field using the arc-length parameterization such 
that g{TT,Tt) = QnvTi^Ti^ = 1, where g is the new metric and tt^ is the momentum 
field redefined in terms of this new parameterization (assuming K-(0) contains no 
derivative couplings). That is, we want to scale our coordinate system in order to 
obtain a conformal transformation of the metric that normalizes our momentum field. 

In fact, there are three possible normalizations, depending on the nature of tt: 



However, this will not concern us here, once our main focus will be to find the 
geodesies of g. Therefore, we have a clear distinction among the three foliations of 
Lorentzian spaces: 

g{'K, tt) > represents the spacelikc-lcaf and contains only non-physical objects — 
after the arc- length reparameterization we will have ^(7f,7r) = 1; 

g{7r,7t) = clearly singular, representing the null-leaf where there is no concept of 
distance, given that all objects here are massless; 




1, spacelike; 
0, hghtlike; 
— 1, timelike. 
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^(•77, tt) < represents the timelike-leaf and contains the physical objects — after 
the arc-length reparameterization we will have g{7r,7t) = —1. 



In order to find the arc-length parameterization, we will borrow an idea from 
classical mechanics called Jacobi's metric [46]. In the Newtonian setting, Jacobi's 
metric gives an intrinsic geometry for the configuration space (resp. phase space), 
where dynamical orbits become geodesies, i.e., it maps every Hamiltonian flow into a 
geodesic one; therefore, solving the equations of motion implies flnding the geodesies 
of Jacobi's metric and vice-versa. This can be done for any closed non-dissipative 
system (with total energy E), regardless of the number of degrees of freedom. 

Before we proceed any further, let us take a look at a couple of simple examples in 
order to motivate our coming definition of Jacobi's metric. 



9.1. Harmonic Oscillator 

Let us consider a model similar to the harmonic oscillator (in (1 + 0)-dimensions): 
without further considerations, we simply allow for the analytic continuation of the 
frequency: /i = +1^"^ or /i = —u^. Its Lagrangian is given by L = ^ =]= ^ g^, where 
fi > 0, from which we conclude that for a fixed E such that E = | (g^ ± /i g^), we have 

So, in the spirit of what was said above, we want to find a reparameterization of the 
time coordinate in order to have the normalization g = 1, where the dot represents a 
derivative with respect to this new time variable. 

Thus, 



dq dq ds 
dt ds dt 



if 



ds 



dt 
.'. ds 



'2[ET^q' 



'2[ET^q' 



dq 
ds 



1; 



'2(ET|g2) dt- 
^ gj, = 2[ET^q']g. 



As expected, the new metric, g^, is a conformal transformation of the original one, 
g; and under this new metric, our original Lagrangian is simply written as. 



9.1. Harmonic Oscillator 
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1 -i -i 2 

= 2 dijQ r T 2 ? ; 

1 A* 2 

= 2£'(^'^) T 2 ^ ' 

Now, as we said above [16], the Euler-Lagrange equations for this conformaUy 
transformed Lagrangian are simply the geodesies of the metric g^. The task of finding 
the geodesies 7(s) of the metric can be more easily accomplished with the help of 
the normalization condition, gsi^j ^) = 1) ^^'^ the initial condition 7(5 = 0) = 0, 
given that, in this fashion, we only need to solve a first order differential equation: 



ds \]2{ETl7') ' 
with the initial condition: 7(0) = . 

Here we should note that solving for ^e(7',7') = 1 (spacelike-leaf) is analogous to 
solving gE{l' = ~1 (timelike-leaf), once the two cases are symmetrical about the 
origin. 

Analytically solving the equation above yields 2 possible answers, depending on the 
particular form of the potential (/i > in both cases): 



1. V+ = +fi 7^/2: For the case of a positive pre-factor, we get that 7 ii{2E — /i7^) '+ 
2 E arcsin (7 ^fJ^2E) - 2 s ^ = 0; and 

2. V- = —jj, 7^/2: For the case of a negative pre-factor, we find that 7 a//U (2 + /i7^) '+ 
2 E arcsinh(7 ^/fi/2E) - 2s ^ = 0. 



It is clear from the expression for V+ that 7^ ^ 2E/fi, i.e., the length of the [classical] 
geodesic is bounded; this does not happen with V_. 

These geodesies, 7-1- , clearly depend on the parameters E and /x; therefore, in order 
to plot 7(s), we have to make 2 distinct choices: E/fi = 1 (left plot) and E/fj, = — 1 
(right plot). The last plot comparatively depicts both geodesies. 
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Figure 1: The top- left-corner plot shows the geodesic for Ej [i= 1, while the top-right-corner one 
depicts it for = —1. The bottom-center plot superimposes both of them. 

This "harmonic oscillator" model already portrays the features that we want to 
identify in forthcoming applications of this technique: changing the values of the 
parameters of the potential we can identify 2 distinct types of geodesies which will be 
related to different solutions on the coming examples. 

Note, however, that in general the parameters will not be able to vary freely: they 
will belong to some specified set; and they will be related to each other, i.e., we will 
be able to find a function of the parameters that constrains their behavior — this will 
establish "dualities" among the parameters of a given theory, as will be shown below. 



9.2. Cubic Potential 



For completeness sakes, let us try another example, that of a cubic potential given 
by: = \ g{q, q) + It's Jacobi Metric is given hy = 2{E + g and the 
equation we need to solve in order to find the geodesies, 7(s), of this metric is. 



gj,ii,i) = 2iE + qy3)gii,i) = l; 

2{E + q'/3) {i{s)Y = l- 
with the initial condition: 7(0) = . 



Once again, we should note that solving for gEil'il') = ~1 (timelike-leaf) is 
analogous to what is being done above: the graphs are reflected with respect to each 
other. 



9.2. Cubic Potential 
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This example is slightly different from the previous one for the following reason: 
before, it was the relative values of E and fi that determined our two solutions, i.e., 
one for E/fi ^ and one for E/^ < 0. Now, this cubic theory has no free parameters 
in its potential, therefore the geodesies are labeled by the arbitrary parameter E. 

Just as before, the metric will be [artificially] degenerate when E = V{'~f), which 
tells us that 7 ^ {—3E)^^^: when E ^ the geodesic is allowed to have any length, 
but when E < the geodesic vanishes for some time, after which it starts to grow. 

This is very interesting because it essentially says that there is a "lag time" before 
this solution comes alive: this leaf is non-existent for some "proper time" and then it 
springs into being quite abruptly. 




Figure 2: The top-left-corner plot shows the geodesic for E ^ 1, while the top-right-corner one 
depicts it for E = —1. The bottom-center plot superimposes both of them. 

Now that we are done with our examples, we are ready to find the generalization of 
Jacobi's metric to the field theoretical setting can be accomplished via the definition 
of the following conformally transformed metric: 



L = ^£?(vr,7r) - V;(0) ; 

where, 

gl = 2{E-VA<P))g; (9.1) 



and E (an arbitrary parameter) is the total energy of the system. 
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Chapter A. \J 

Applications in Quantum Field Theory 



We will use the approach in terms of Feynman Path Integrals in order to make things 
more straightforward. However, we could as well talk in terms of the phase space of a 
given QFT and its vacuum manifold, i.e., its moduli space. 

Starting from the partition function, we have the following (h — l): 

Z[J] = Je^5M+^/-^W'^W'^''^D0 ; 

= X J" e* f 9Ei^,Tr) + J(x) <t>i^) d'^x rj^^ . 

where K is a normalization constant such that Z[J = 0] = 1, and we have already 
written the Action in terms of Jacobi's metric. 

Now, let us expand the partition function above in terms of its classical part and 
its quantum fluctuations, i.e., = 0ci + ^4>, and the classical Action is given by 
Sci = S[(f)c\] = gEi'^cij'^ci), where M is the particular region of spacetime where 
the integration is performed and tt — d(f) — for this calculation we are assuming the 
fluctuations vanish at the boundary, 50|afw = 0. Thus, we have that. 



Z = Jexp{i5[0ci + 50]}2)0; 

= Jexp|i j ^E(7rci,7rci) +^E((^7f,(57f)| ; (10.1) 
^e'^'^'FridM]; (10.2) 

where, 

Fr[dM] = r e^^l'^'^] D((50) = r f^a^^is^'^^) D((50) ~ ^ ■ (10.3) 

Vdet(-D2) 

where is the [quantum] fluctuation part of the partition function, the index 

T denoting collectively the parameters of the potential (mass, coupling constants. 
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10. Applications in Quantum Field Theory 



etc), and = is the Laplace-Beltrami operator constructed from the covari- 

ant derivative (V^) associated to the Levi-Civita connection of There are two 
important things to note from (10.3): the fluctuation term can only depend on the 
parameters of the potential and on the boundary dM, i.e., surface terms (that we 
assumed vanishing); and its structure is analogous to that of the classical part, in that 
is the same on both. 

However, the above presumes that there is only one solution to the theory in question. 
But, the method we are developing is exactly to use Jacobi's metric in order to find 
and classify all of the solutions of the theory given, i.e., we want to be able to use 
this tool to study the moduli space of the problem at hand. Therefore, we need to 
generalize the situation above for the case of many solutions, which is not a difficult 
task: 



z=n 



TV 



e'^^iF^[dM]; (10.4) 

-1—1- p'-^cl 



det{-Dl) 



where u counts the number of different solutions (i.e., the number of critical points 
of S'[0], including its multiplicity) denoted by (p'^i, and S^i = S[(f)'^i]. Note that, in 
order to find all of the possible critical points of S[(l)], we need to take into account its 
r-dependence, and in doing so we are implicitly assuming that r G C, ensuring we are 
able to find them all. 

There are two major observations to be done at this point: 

1. There are two kinds of discontinuities present in the above construction, [ ]: first, 
critical points can coalesce, which happens in the complexified catastrophe set 
in parameter space; second, critical values 53(iS'[0J^J) = S^i can coalesce, which 
happens on the Stokes set in parameter space and corresponds to the appearance 
or disappearance of a subdominant exponential in a "non-local bifurcation". 
These correspond, respectively, to realizing that the partition function is to 
be taken over complex fields (rather than real ones, as is customary), where 
the contour of integration (rendering the partition function finite, [30]) will 
ultimately determine the parameter space (these are the Lee- Yang zeros of 
our theory); and the coalescing of S^y represents the Stokes phenomena of our 
theory. Both of these will determine the phase structure of our problem (see 
[35]). Note that it is the contour of integration that connects both of these, once 
the appropriately chosen range of will render e^'^t'^^ convergent. 

2. The asymptotic behavior of the partition function depends only on the critical 
points of S, i.e., on (p'^y. Thus, when 5'[0] (resp. ^[0]) is a Morse function, i.e, 
smooth with no degenerate critical points, we can use Morse's lemma to show 
that the critical points are isolated, keeping in mind that the number of isolated 
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points is a topological invariant. (It is worth noting this can be generalized in 
the thermodynamical limit via the Morse-Palais lemma, just as we can relax 
the condition of non-degeneracy of the critical points via Morse-Bott theory.) 
Using this, it can be shown that 0-space (resp. Phase Space) is a CW-complex 
with a i>-cel\ for each critical point of index u: the fluctuation term contains the 
Maslov-Morse index (in by means of its dependence on K-[(^]) that accounts 
for the discretization of Path Space and corrects for the thermodynamic limit 
of the particular solution in question, i.e., every time the denominator in (10.5) 
vanishes, F^[d9i{] passes a singularity in such a way as to ensure the proper 
phase; the phase factor that arises in this way is nothing but e*"*^, where a is 
some angle and u counts the number of zeros (with multiplicity) encountered 
along the particular path in question — it is called the Maslov-Morse index. 

Now, we are ready to employ this machinery in order to draw several important 
conclusions. Here are them: 

1. The partition function, seen as a function of the parameters of the potential, Z = 
Z[t] = 2.[mass, coupling constants], is a meromorphic function: it is holomorphic 
on a subset of C except for the set of isolated points given by the values of the 
parameters along Stokes' lines (resp. critical lines of phase transition). We have 
to keep in mind that the parameter space had to be complexified in order to 
yield all possible solutions to the theory given; another way to think about this 
is in terms of the contours that render the partition function finite: changing 
these contours (looking for all possible ones that make the partition function 
converge) will affect the allowed values for the parameters, i.e., these contours 
ultimately determine the parameter space, as done in [.■)■">] (see also [10]). 

2. Under the [full] elliptic modular group, F = SL(2,Z) = {{cd) \ a = ^, b = 
0, c = 0, d = ij (see [38]), the partition function is a modular function, i.e., 
for any M e F we have that 2, [Mr] = Z[t], where Mr = f^. (For a fuller 
appreciation of the importance of this fact, see, e.g., [39].) 

3. The Action, written in terms of Jacobi's metric, S[(f)] = J^gEi'^^^)^ can be 
thought of as the Morse-theoretic Energy functional, £'[7] = / g{'~^','j'), for 
the path 7, where 7' = d'j/ds. Therefore, we can readily say that all of the 
critical points of S[(j)] are given by minimal [Lagrangian] manifolds, and, just like 
E", S" = 5(j){x)5(j){y) is a well defined symmetric bilinear functional. This 
implies that S" = if, and only if, vr is a Jacobi field, which, in turn, implies that 
^IdM = 0. Therefore, using the split of the partition function (10.4) in terms of 
a classical part and its quantum fiuctuations, we see that the vacuum manifold 
(moduli space) of the quantum theory is given by its classical minimal Lagrangian 
manifold with the quantum corrections being given by extensions of it via Jacobi 
fields. That is, the classical minimal manifold is extended via the gluing of the 
quantum fluctuation manifold obtained as a solution to the equations of motion 
coming from Sa = /^^^(^vr, (57r). Thus, the quantum corrections are handles 
attached to the classical solution, without changing the classical topology. 
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4. Continuing the reasoning above, it is not difficult to see that a phase transition 
will happen when the quantum corrections (in terms of handle attachments) 
actually do change the topology of the particular solution in question. In fact, 
if V^[0] crosses a critical point of index z/, the handle to be attached is a u-cell 
(i.e., a i/-simplex) — the connection with what has been said before should be 
straightforward. Analogously, we can say the following: let = V~^[E] = 
{(j) I = E}, i.e., Tif is the set of all field configurations which have the same 
potential energy — this manifold is the configuration space (resp. moduli space), 
an equipotential surface. Thus, the family of equipotentials {T,^}Em foliates the 
configuration space (moduli space) in such a way that if S;^ is '^°°-diffeomorphic 
to Ef (for two different values of the energy, E and E) then there is no phase 
transition. Conversely, a phase transition will be characterized by the existence 
of a certain critical value Ec such that {'E^}e<Ec is not ^°°-diffeomorphic to 
{Ef }e>£;^ (note that, upon a phase transition, the parameters of the potential 
change from r to r). Further, the topological difference between these two leaves 
is a u-cell, where u is the index of the critical point E^. Loosely speaking, it can 
be said that the origin of phase transitions is this change in topology. Therefore, 
different phases of the theory are topologically inequivalent. 

5. Finally, let us note that all of these conclusions we have drawn so far have one 
last implication, a topological constraint: the Euler Characteristic, x, computed 
from gives the quantization rules [for the energy]. Therefore, our theory can 
be topologically quantized; in fact, each topologically inequivalent leaf {Ef } has 
its own Xt ^^d, thus, its own quantization rules. 

With all of these facts in hand, let us see what is the framework they imply: 
given a theory, L = gt(7r,7r) + Vr[4>], we can readily construct its Jacobi metric, 
g'^ = 2 {E — Vr) g: from this, we can do two things, either compute the geodesies of 
and classify their Jacobi fields in terms of r, or calculate its Euler Characteristics, 
Xri aiid find the quantization rules (bearing in mind that they will vary with r). The 
different geodesies (corresponding to different Jacobi fields) will label diffeomorphically 
equivalent foliations {Ef }, while the order of the zeros of F^[dM] will correspond to 
the z/-cells associated with a particular phase transition — furthermore, they will be 
responsible for the Lee- Yang zeros and the Stokes phenomena of the theory. In turn, 
this gives the partition function a meromorphic and a modular character when seen 
with respect to r. 

All of this was possible because we extended the parameter space, allowing r to be 
complex. This is completely analogous to complexifying the solution space (moduli 
space) of the theory given, a fact which is made clear in [35] (and employed in lattice 
calculations in [ ' ]). 

Lastly, using r, it will be possible to construct "dualities" between different phases. 
However, in general, these dualities will be non-trivial combinations of the parameters 
of the theory (as opposed to what we found in the simple examples above, where 
E/ fi —E/ji or E —E did the job). 



10.1. The A 0* Potential 



51 



10.1. The Potential 

This theory is defined for scalar-valued fields, 0, by the Lagrangian L — ^ {g{T^-iT^) — 
jicf)^ — ^ (p^, i.e., the potential is given by — | (//^^ + | 0"^); which is invariant 
by Z2-reflection: i— > —(f). 

The Jacobi metric for this theory is given by = 2 (i? — ^ 0^ — | 0"^) g. And 
solving the normalization condition g'^il',!') = 1 — where 7(5) is the geodesic we 
want to compute and 7' = where s is the arc- length parameter — should give us 
the geodesic structure of the vacuum manifold (resp. moduli space): 

^^(7,7) = 1; 
2(£;-ff -t7')(7')' = i; (10.6) 

The full picture presents itself upon a more detailed analysis of j {E ~ /U7^/2 — 

1 /2 

A 7^/4) (^7, which is the [elliptic] integral that needs to be solved in order to find 
7(s). Therefore, it is useful to consider the polynomial -P(7) = (7^ — ri) (7^ — r2), 
where (-A/4) P(7) = E - ii^^ jl - A 774 and ri,2 = -(/x ± ^^EX^^i^ ) /A. It is 
the discriminant of this polynomial -P(7) that will, ultimately, determine the different 
solutions of the theory: A = (ri — r2)^ = XE + iJ?/A ^ 0. 

As mentioned above, the "dualities" between different phases involves non-trivial 
combinations of the parameters of the theory. In this case, the dualities are given by 
the possible values of the discriminant above: A > 0, A = or A < 0. 

The two inequalities, A > and A < 0, are related by the reflection of the r 
parameter of this theory, where r = /^^/A, i.e., by the analytic continuation of the 
mass parameter such that //^ 1— > — //^, which implies that r 1— > — r — note that 
this can be obtained by a modular transformation, Mr = = — r, such that 

6 = = c, (i=l and a = —1; further, the analytic continuation of /j, is such that 
fj, ^ ±i fi, or more gcncrically as 1— > e^*'^'^/^ /x, for odd v (this is related to the Stokes 
phenomena discussed previously, where 1/ selects a certain Riemann sheet for this 
analytic continuation, which depends on the Maslov-Morsc index u) — ; meanwhile, the 
case A = has to be computed separately, once it establishes a fixed value r = —4 E 
— note that when A = the resultant of P{'~f) and its derivative, P'i'j), also vanishes, 
this resultant being given by Res(P, P') = — | A, which says that P and P' have a 

2 

common root (which happens when t = ^ = —4:E). 

This whole framework is potentially a very interesting result, once discriminants 
arc closely related and contain information about ramifications ( "branching out" or 
"branches coming together") in Number Theory. As we just saw, these ramifications 
are related to the topological structure of each solution, where the Maslov-Morse 
index, u, selects the appropriate Riemann sheet (for the analytic continuation) and 
the topology of the problem (determined by the attachment of a u-celV). Therefore, 
this duality relating the two different solutions, for r and — r, is given by a [particular] 
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modular transformation of r (resp. analytic continuation of fi), which turns out to be 
a measure of the ramifications of the problem. 

The analytical answers for the geodesic 7(5) are found to be given by the following 
implicit equations: 



A = 0, Res(P,P') = : 



1 o u 



A > 0, Res(P,P') < : 

3 s + 7 a/ (7^ - ri) (72 - + 2 ri y/r^ F 



7 



2 fi 

3 A 



7 



E 



7 



0; 



(10.7) 



(10.^ 



A < 0, Res(P,P') > : 

3 s + 7 a/(72 - ri) (72 - r2) + 2 ri ^/r^ F ^ 



2 fi 

3 A 



7 



7 



10.9) 



. 



Note that, in the equations above, F{z; k) is the incomplete elliptic integral of the 
first kind, while E{z; k) is the incomplete elliptic integral of the second kind. Moreover, 
ri^2 are defined (as shown above) as the solutions to the polynomial equation ^(7) = 0: 
ri2 = -(/i ± + 4:XE)/X. On top of this, although (10.8) and (10.9) have the 
same form, they will yield distinct solutions, once the relation given by A = XE + fi'^/4 
will either be positive or negative, which, in turn, affects the outcome of ri 2 — as 
already discussed above for r — r. 

The graphical results are shown below, where y = 1,7(0) = and A, respectively, 
assumes positive (A > 0), null (A = 0) and negative (A < 0) values: 
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We can clearly see that, for positive A (leftmost graph), we have a smooth geodesic 
7(s) representing the symmetric phase. Then, when A vanishes (middle graph), there 
is a clear change which delimits the 2 different phases of the theory. Moreover, when 
A is negative (rightmost graph) , we have the third [broken-symmetric] phase of the 
theory (which has a finite geodesic). 
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Chapter A. A 

Gauge Theory Examples 



Now, we will consider examples of theories with gauge symmetry (rather than the 
discrete Z2-symmetry of above). As seen in (10.6), taking 7 1— >■ —7 did not change the 
geodesic structure, i.e., the procedure was covariant (resp. equivariant). This should 
come as no surprise, once the potential, K-(0), is an equivariant function of the fields 
and the spacetime (base manifold) metric does not depend on the gauge symmetry 
involved either. 

As mentioned before in section 10, we have assumed thus far that the quantum 
fluctuations vanish at the boundary, 50|a5w = 0. This, in turn, implies that the Jacobi 
fields vr will satisfy tcIs^ = 0. While the fields considered have no internal structure 
(i.e., gauge symmetry), this is a fairly straightforward constraint, in the sense that its 
solution is trivial (albeit labeled by r). However, when gauge symmetry is present, 
there may be non-trivial solutions to these constraints. This is the phenomenon of 
spontaneous symmetry breaking. 

Therefore, following the discussion done in section 10, the following may happen: the 
quantum corrections, which are handles attached to the classical solution, may now be 
able to change the topology of the classical solution. So, when a symmetry is reduced 
from G H (where if is a subgroup of G; see, e.g., [25, 26, 27, 28, 29, 31]), we start 
with a Jacobi field which satisfies ttgIssw = and end up with a Jacobi field who 
only has H as symmetry, which implies that TiH\d9^ = — this means that the initial 
degrees-of-freedom that combined (respecting the G symmetry) to yield the [initial] 
constraint TiG\d(M = 0, are not all available now, such that only part of the original 
symmetry is still respected, yielding TTHlas^ = (the remaining degrees-of-freedom 
having recombined in non-trivial ways); thus, we end up with a source (or sink) of 
Jacobi fields that only have H as symmetry. 

In this sense, the Atiyah-Singer Index Theorem can be used to measure this variation 
(i.e., to measure the inequivalent representations of the algebra of observables): 

1. If all of the quantum corrections preserve the topology of the classical solution, 
the topological index does not change, which implies that the analytical index of 
the differential operator in question also does not change, which, in turn, leaves 
the vacuum state unchanged. 
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2. On the other hand, if the quantum fluctuations change the topology of the 
classical solution (as described above), the topological index will change (following 
the attachment of the appropriate z/-cell), implying that the analytical index 
changes as well, which means that the vacuum state changes. 

From a different viewpoint, the question can be posed in the following way: Given 
a certain symmetry breaking connection, how can the topology of the moduli space of 
the associated Higgs Bundle be studied? 

There are some studies in this direction, but no general answers: this is because 
these types of characterizations are highly model-dependent, i.e., they depend on the 
particular properties of the connection chosen for the Higgs Bundle in question. For 
instance, following the discussion above, if the boundary DM is non-existent (i.e., M 
is compact), then 7f|a5w = is trivialy satisfied; however, if M has, e.g., punctures 
(i.e., a finite set of points omitted), then the connection in question will have a certain 
ramification structure and the constraint vf |a5w = will have non-trivial solutions (see, 
e.g., [52], and references therein); therefore, the boundary conditions (in the sense of 
[35, 40]) on the Jacobi fields ultimately determine the structure of the connection (i.e., 
gauge field). 

The examples below are twofold: the first one ( Landau- Ginzburg Functional, section 
11.1) serves the purpose of showing the equivariance of the method discussed in this 
work, once the connection (which is u(l)-valued) does not undergo symmetry breaking; 
while the second one (Seiberg-Witten Functional, section 11.2) can be understood 
as defined for pairs (A, 0), where A is a Hermitian connection (compatible with the 
holomorphic structure on the bundle in question) and is a section (of the bundle at 
hand) — in this sense, one can consider the space of solutions to Hitchin's equations 
(i.e., the moduli space of the associated Higgs Bundle) given by. 



where Fa is the curvature of A and d'^ip is the anti-holomorphic part of the covariant 
derivative of (compare these with equations (11.11) and (11.12)). In this sense, the 
different solutions found via (11.13) are a direct statement about the topology of the 
moduli space of the Higgs Bundle under study. 



To start off, let us consider the case where the base manifold is a compact Riemann 
surface S equipped with a conformal metric and the vector bundle is a Hermitian line 
bundle L (i.e., with fiber C and a Hermitian metric (■, ■) on the fibers). 

The Landau- Ginzburg functional is defined for a section and a unitary connection 
Da = d+y4 of X as ((T G M is a real scalar). 



F^ + [0,0*] = O; 
c/'^0 = ; 



(11.1) 
(11.2) 



11.1. The Landau-Ginzburg Functional 



11.1. The Landau-Ginzburg Functional 
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Thus, its Euler-Lagrange equations [of motion] are given by: 

D\Fa = -Re{DAif,if) ; 

where, D\ is the dual of Da, i.e., D\ = — * Da* = — * (d+A) *. Note that the 
equation for Fa (the second one above) is Unear in A. Since Da is a unitary connection, 
A is a u(l)-valued 1-form. This Lie algebra (of the group U{1)) will sometimes be 
identified with iM. — in other words, our A corresponds to —i A in the standard 
physics literature (where A is real- valued) . 

Before we go any further, some notational remarks are in order. We decompose the 
space of 1-forms, on E as = fi^'*^©!^"'^, with spanned by 1-forms of the type 
dz and by 1-forms of the type dz. Here z = x + iy is a. local conformal parameter 
on E and z = x — iy. Therefore, dz = dx + i dy, dz = dx — i dy, = \ {d^ — i dy) and 
dz = \ {dx + i dy). Furthermore, if dx and dy are an orthonormal basis of the tangent 
space of E at a given point, we have that {dz, dz) = 2, {dz, dz) = 2 and {dz, dz) = 0. 
Given that the decomposition = Q^'^ © Q^'^ is orthogonal, we may also decompose 
Da accordingly: Da = 9a + Ba, where dA^ G fi^'°(x) and Ba^ G fi°'^(x) for all 
sections of X (holomorphic, dAf{z,z) = -v^ f{z,z) = f{z), and anti-holomorphic, 
dAf{z,z) = -v^ f{.z,z) = f{z), parts; i.e., the space of 1-forms and the space of 
connections is decomposable into a direct sum of its holomorphic and anti-holomorphic 
parts: dA = d + and Ba = d + A^'^; while the exterior derivative is given by 
(i = d + B). As expected, we have that 9^ <9a = = 9^ Ba and Fa = — (c^a + ^^a) • 

It is not difficult to show [HO] that: 



L{^,A)= I \FA? + \DAV? + \{a-M''f ; 

= 27r deg(i:) + Jj\Ba^\^ + {*HFa) - I ~ 1^1'))' ; 

1 

where deg(x) = ci{l) = — tr(FA) ; 

Z TT 

i.e., the degree of the line bundle is given by the 1st Chern class. 

Therefore, a very useful consequence from the above is that if deg(x) ^ the lowest 
possible value for L{ip, A) is realized if and A satisfy 



BaV 



0; 
1 



*(^FA) = -(a-|y.|2) . 



:ii.3) 
:ii.4) 
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These equations are just the expression of the self-duahty of the Landau- Ginzburg 
functional. If deg(x) < 0, then these equations cannot have any solution, thus one has 
to consider the self-duality equations arising from the Landau-Ginzburg functional 
where the term +2 7r deg(x) is substituted by —2-n deg(x). Therefore, without loss 
of generality, we shall assume deg(x) ^ 0. A necessary condition for the solvability of 
(11.4) is that. 



2 TT deg(x) = j iFa = \ j {<J- Iv^p) ^ I Area(S) ; 
. 47rdeg(x) 

and the equahty only occurs if, and only if, = 0. 

The following result is useful when studying the solutions of the above functionals 
[50]: Let be a compact Riemann surface with a conformal metric and L as before. For 
any solution of (11.3), we have that \(p\ ^ a on S. That is, this maximum principle 
states that the amplitude of the field cannot exceed the height of the potential (see 
the first plot below). 

Let us now construct the Jacobi metric for this potential and find its possible 
geodesies: 



gE = 2{E-V,{j))g; 



In a complete analogy to what was previously done, we consider the -P(7) = E — V{'~f) 
polynomial. However, in order to make this analysis clearer, let us, in fact, use a 
slightly different polynomial given by -P'(7) = 4 (-P(7) — E^. It is straightforward to 
see that P"(7) — (|7P~'^) ~ ^')^ which means that a is the only root of P'{j), 
with double multiplicity (the two roots coalesce into one). 

This situation implies that the discriminant of P'{'y) vanishes, i.e., A = 0, and the 
different phases of the theory are labelled by a = and a > 0. 

The plots below have, respectively, the following values for a: 0.0, 3.0, 5.0, 7.0, 11.0 
and 31.0. As they show, when a = we have Ci(x) = deg(x) = 2^ tr(F^) = and 
its geodesic has a clear character which is quite different otherwise: 
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Figure 4: The plots above show the clear distinction between the a — solution and the a > 
ones; the first plot, on the upper left corner, has cr = 0, while the other ones have differing positive 

values for a. 

As a last remark on this model, note it can be thought as defined by a pair [A, 
where A is a unitary connection compatible with the holomorphic structure of the 
bundle in question, and is a [global] section of this bundle. In this sense, Hitchin's 
self-duality equations become, 

dA^ = 0; (11.6) 
Fa+[^,^*]=0. (11.7) 

Compare these with (11.3) and (11.4). 

In this sense, the solutions found above are a direct statement about the topology 
of the Higgs Bundle in question. 



11.2. The Seiberg-Witten Functional 

In this case, the base manifold, is a compact, oriented, 4-dimensional Riemannian 
manifold endowed with a spin^ structure, i.e., a spin'^ manifold. The determinant line 
of this spin*^ structure will be denoted by L and the Dirac operator determined by a 
unitary connection A on L will be denoted by "Da- Recalling the half spin bundle §^ 
defined by the spin'^ structure, we see that 2)^ maps sections of into sections of S''' 
[47, 48, 49, 50]. 

In this fashion the Seiberg-Witten functional for a unitary connection A on L and 
a section ip of S"^ is given by, 

SW[^, A] = \VaV? + WX? + + (11.8) 

where Va is the spin^ connection induced by A and the Levi-Civita connection of 

is the self-dual part of the curvature of A and R is the scalar curvature of M. Its 
Euler-Lagrange equations are given by. 



V*A Va^ 

d* Ft 



R 1 , , 
4+4"^" 



(11.9) 
11.10) 
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Using a spin frame, it is not difficult 
can be written in the following form: 



to show that the Seiberg-Witten functional 



SW[ip,A] 



5W 



^1, J, k 



A. As a 



where are 1-forms dual to the tangent vectors e^, e^{ek) 
corollary of the above, the lowest possible value of the Seiberg-Witten functional is 
achieved if ip and A are solutions of the Seiberg-Witten equations: 



F 



:il.ll) 
:il.l2) 



Thus, self-duality is at work yet again: the absolute minima of the Seiberg-Witten 
functional satisfy not only the second order equations (11.9) and (11.10), but also the 
first order Seiberg-Witten equations (11.11) and (11.12). 

Although our discussion of the Seiberg-Witten functional, so far, has mirrored our 
discussion of the Landau- Ginz burg one, the parameter a on the latter has had no 
analogue in the former. This can be accomplished with the introduction of a 2-form /i 
and the consideration of the perturbed functional, 



"a - ^ (^i ■ Cfc ■ V9, v?) e-' A e + /i 



/ \VAP\' + \Fi\' + 



R 



+ 2(F+,/i) 



Their corresponding first order equations of motion are. 



VaP = 

=^ (cj ■ Ck ■ p,p) Ae'' - fi . 

If we assume that is closed and self-dual, then we see that {FA,fi) = {F^,fi), 
once {F^,fj,) = due to the orthogonality between anti-self-dual and self-dual forms. 
Thus, given that Fa represents the first Chern class ci(x) of the line bundle L, and 
we assumed /i to be closed (hence it represents a cohomology class [/i]), the integral 

does not depend on the connection A, thus representing a topological invariant, denoted 
by {c^{L)A[^J])[^^{]. 

Just as before, we also have a maximum principle: For any solution p of (11.9) — 
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in particular, for any solution of (11.11) 
manifold, we have that. 



on a compact ^-dimensional Riemannian 



max |</)p ^ max( 



-i?(x),0) . 



As a direct consequence of this, if the compact, oriented, Riemannian Spin'^ manifold 
M has nonnegative scalar curvature, the only possible solution of the Seiberg-Witten 
equations is. 



^ = ; F+ = . 

The Jacobi metric for this Seiberg-Witten model is given by 

1 , 



9e 



171 



9 ; 



(11.13) 

for a geodesic 7; and, just like before, although ^(7) = E—V{^) = E+R |7p/4+|7|''/8, 
let us consider P'(7) = 8 (P(7) - E) = (I7P - O) (I7P + 2R). It is then clear that 
when R = the discriminant vanishes (A = 0) and the 2 roots merge into 1, what 
constitutes one of the phases of the theory. When R ^ 0, the discriminant is either 
A > (-R > 0) or A < (i? < 0), which accounts for the other 2 phases of the theory. 

The plots below were obtained with the choice of i^^ = 1.0 and R respectively equal 
to 0.0, 3.0, 7.0, -3.0, -5.0, -7.0. 



1 






i 




L — ^"^1 


r 

.... _.L/ 


! 




/\ 


1 


1 


/ 


1 


1 










3£ = S[E4-«Fin(= + |yl'ja)B 































T 1 














- - ■- 


_ j_ _i _ _ _ __ 


























-/- 


1 





















/ 










_ 






















- y - — - 
















/ 







Figure 5: The plot in the upper left corner depicts the solution for i? = 0, while the plots in the 
upper center and upper right corner show solutions for i? > 0. The plots in the lower row picture the 

solutions for i? < 0. 
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Chapter 

Conclusions 



Using a conformal transformation that amounts to finding the arc-length reparameter- 
ization of the given problem, we were able to construct a new metric, called Jacobi 
metric, such that its geodesic equation is equivalent to the original equations of motion. 

Then, we separated the given problem into its classical part and quantum fluctuations 
and, by realizing that the action written in terms of the Jacobi metric is the so-called 
"energy" in Morse Theory (resp. Morse-Bott theory), we see can reinterpret this as 
a "topological expansion" , in the sense that the quantum fluctuations are handles 
attached to the classical solution. 

Further, when these handles change the original (classical) topology, via the gluing 
of an appropriate z/-cell, there is a phase transition, in the sense that we move from 
one solution [of our QFT] to another. This can be measured via the use of an Index 
Theorem: once it relates the analytical index of the differential operator in question 
with the topological index of the manifold under study, we clearly see that if the 
topological index changes (because of a certain handle attachment), the analytical 
index must change as well; which means that the zero modes, the solutions of the 
equations of motion, have to change. 

Therefore, by studying the geodesies obtained from the Jacobi metric with respect 
to its parameters, we can compute the topological index (via generalizations of the 
Gauss-Bonnet theorem, using the Jacobi metric we derived before) and thus compute 
the zero-mode solutions with respect to the different values of the parameters. In this 
way, we are able to classify all possible solutions to the QFT in question, keeping in 
mind that the quantum corrections are given by handle attachments that may or may 
not (depending on the values of the parameters) change the topology of the classical 
solutions. 

To this picture, we add Lee- Yang zeros and Stokes phenomena, in order to obtain a 
more robust view of what is at stake: the Lee- Yang zeros accumulate along Stokes 
lines and pinch the parameter space, creating different regions of "allowed values" 
for the parameters of the theory (mass, coupling constants, etc). This represents 
different sectors of theory, different phases of the theory. And we showed that these 
are topologically inequivalent. 

In this fashion, the partition function is meromorphic, once it is singular along the 
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12. Conclusions 



Stokes' lines. However, it has various degrees of modular symmetry, depending on 
the particulars of the theory in question. And, in this sense, different sectors of the 
theory are related to each other [by an appropriate modular transformation]. 

In turn, this implies that different solutions of the theory are related to each other, 
a fact that we dubbed "duality" . Ultimately, these dualities are determined by the 
actual values of the parameters (and how they compose in order to create the modular 
symmetry in question), which are determined by the boundary conditions of the 
Schwinger-Dyson equations (as explained earlier). 

Future work will focus on D-modules and dimensional construction of 0-dimensional 
theories (which can be completely solved). We also intend to generalize the cubic 
potential showed earlier to matrix- and Lie-algebra- valued fields [j j], analytically 
solving the its 0-dimensional counterpart and dimensionally constructing it (via D- 
modules): it seems plausible that such an extension might be related to current 
developments in three-dimensional gravity, [ ]. 



Part III. 

Three-Dimensional Gravity and Airy 

Functions 
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Chapter 



13 



Introduction 



Recently, new developments have shone new light on the problem of three (1 temporal, 
2 spatial) dimensional gravity, [39]. Loosely speaking, here is the gist of the matter. 

Classically, 3-dimensional gravity (plus a cosmological constant) is given by the 
action. 



It is not difficult to see that the solutions to the above are all locally equivalent and 
there are no gravitational waves. Thus, naively, one may think that this problem is 
tractable quantum mechanically. 

However, on second thought, power counting tells us that G has dimensions of 
length, rendering the theory unrenormalizable. This leads to the conclusion that the 
quantum theory does not exist. But this is premature as well, once divergences in 
perturbation theory can be removed by field redefinitions {g^j^y ^ g^u + a R^j_y + ■ ■ ■) 
and a renormalization of 

Therefore, a slightly different venue of attack must be drawn in order to resolve the 
conflicts above. What is customarily done is to write the above Action as a gauge 
theory, which will lead us towards a Chern-Simons theory, in the following way: write 
the metric in terms of its Cartan connection, gfiu = Vp-i^ ^^/i ^ii known as the "vielbein" 
(triad). From this, we write the spin connection, which is an so(2, l)-valued 1-form: 



At this point, we are ready to combine the vielbein and the spin connection into an 
so(2, 2)-valued connection. A: 



The Chern-Simons action, which is the gauge theoretical version of this problem, is 
thus given by. 




,a b 






67 



68 
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Now, we can try and tackle this problem perturbatively. What we see is that, once 
there is no appearance of the metric tensor in the action above, there are no local 
densities in this problem, so there are no counterterm contributions in this sense. 

Further, we realize that the cosmological constant plays the role of the structure 
constant of the gauge group. Therefore, it is not renormalizable. 

This Chern-Simons description is valid only if the vielbein is invertible, which is 
true for classical solutions; perturbation theory will not change this. 

Finally, this theory is finite and renormalizable by power counting, once there are 
no local counterterms. Therefore, the quantum theory seems to exist. 

However, let us make some non-perturbative considerations: 

• There is some unclear relationship between non-invertible vielbeins and the 
above picture for 3-dim gravity: classical solutions analogous to A — u; — e — 0, 
which are clearly non-geometrical, must be included in a quantized theory of 
3-dim gravity. 

• The equivalence between diffeomorphisms and gauge transformations is not 
straightforward, in the sense that gauge transformations arc continuously con- 
nected to the identity, while more general diffeomorphisms also play a role in 
3-dim gravity. There is no "natural" way to address this situation. 

• Further, the above picture in terms of a Chcrn-Simons theory does not require a 
sum over topologies, a basic ingredient of a theory of quantum gravity. There is 
no a priori reason to add this to the picture described so far. 

• Lastly, for negative values of the cosmological constant, there is a whole class of 
classical solutions known as BTZ black holes (which make their appearance in 
the AdS/CFT correspondence). This clearly shows that the theory cannot be 
"trivial" . 

Therefore, we can summarize the above results in terms of the cosmological constant, 
A, in the following fashion: 

A > : it is not possible to define any set of observables that can be intrinsically 
measured. 

A = : this theory has no gravitons and no black holes, thus it has no S-matrix (no 
degrees-of- freedom) . 

A < : this theory contains BTZ black holes as degrees-of-freedom; therefore its dual 
2-dim CFT can be computed via the AdS/CFT duality. 

So, at this point, the plan shifts to being that of finding the 2-dim CFT dual (via 
AdS/CFT) to the Chern-Simons theory above, with negative cosmological constant; 
for we know that there are degrees-of-freedom in this theory, given by the BTZ black 
holes. 

In this sense, the Euclidean partition function can be written as. 
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Z[I3, e] = tr exp{-p H-iOJ}; 

where H is the Hamiltonian, J is the angular momentum (rotation of the asymptotic 
AdSa), and (3 is the imaginary time. This integral is taken over 3-geometries conformal 
at infinity to a 2-torus with modular parameter t = 6/27i + i f3. 
However, there are two problems with this picture: 

1. The sum of known contributions is not physically sensible because it cannot be 
written as tr exp{— /3 H — i9 J} such that [H, J] = 0. 

2. This partition function is not convergent, once the action is not bounded from 
below. 

To remedy these, we note the following [39]: 

• Real saddle-points do not account for the whole theory. 

• Complex saddle-points are needed, which lead us towards a holomorphic factor- 
ized partition function consistent with an interpretation as tr exp{— /? H — iO J}. 

• Classical geometry: in the semiclassical limit, where G ^ with fixed AdS 
radius we obtain complex (3 and — justifying the complex saddle-points 
above. 

• Finally, the Hawking-Page phase transition with respect to /3, between a thermal 
gas and a black hole, amounts to being a condensation, on the phase boundary, 
of the Lee- Yang zeros of the partition function. 

It is at this point that we make contact with what was said in Part II of this 
thesis, and also in [12, 13, 40, 41]: the parameter space (i.e., the (/3, ^)-space) must be 
complexified in order to account for all possible solutions of the theory — this implies 
that complex geometries must be considered in order to make sense of this theory. 

Furthermore, all ingredients already described in Part II make their appearance 
here as well: Lee- Yang zeros accumulate along Stokes' lines in order to account for 
the Hawking-Page phase transition mentioned, and this corresponds to a topological 
transition between the two phases. Also, in this case, the partition function has 
modular symmetry, Z\t\ = Z[M.t], where M G SL(2,C)/Z2. 

So, let us see what we can say about this in the next chapters. 
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Extended Airy Functions 



The first thing we notice is that the Airy functions (Ai and Bi) can be extended from 
scalar- valued to matrix- and lie-algebra- valued functions, as done in [5:')]. 
In this sense, we end up with an integral representation of the form. 



where V can be either the space of n x n Hermitian matrices (for the case of a matrix- 
valued Airy function) or, in general, a vector space that contains the representation of 
the Lie algebra in question (for the case of a Lie-algebra- valued Airy function); and 



This integral does not exist in the sense of Lebesgue, but it is well-defined as a 
distribution. The formal similarity of these extensions of the Airy function to the 
scalar-valued Airy function accounts for its name. 

In this sense, just as remarked before in chapter 9.2 of Part II, an integral that, 
a priori, looks ill-defined (for it is unbounded from below), does indeed have a 
finite solution and non-trivial content, given by the Airy function (and its respective 
extensions). However, this extensions work for more general polynomials of as long 
as these have the so-called Airy property, ['>:)]. In particular, the differential equations 
associated with this kind of integral representations is treated by Harish-Chandra's 
study of invariant differential operators and integrals on semisimple Lie algebras, [54]. 
Unfortunately, the study of the connections between Harish-Chandra theory, [ i J], and 
the Schwinger-Dyson equations of these models is beyond the scope of this work. 

It is worth remarking that these extensions of $ and J are completely analogous to 
what we did before when we complexified the 0-space. In fact, this is a generalization 
of that idea: you can extend from M-valued to C-valued scalars, but also to more 
complicated fields, as matrix- or Lie-algebra-valued ones. 

Once again, as observed before, the allowed values of J are determined by the 
particular contours that render the integral representation above finite, and there will 
be a certain "duality" between the different solutions determined by the distinct set of 
allowed values (which is, in turn, determined by the particular contour in question). 

Therefore, summarizing, we have managed to generalize the original cubic potential 





J G V. 
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in different ways: from M- valued to C- valued scalars, and from scalars to matrix- and 
Lie-algebra-valued fields, via equation (14.1). Now, two points are important to note: 

1. Equation (14.1) is just a particular case of a more general one, given by [53]: 



where p($) is a G-invariant polynomial on V, and G is the Lie group in question. 
Thus, the Airy case is recovered when we specialize to the case of p{^) = 
$^/3; but other generalizations are possible, e.g., p($) = tr$^/4, or p($) = 
tr (A $ V2 + B $ V4) , where A, B G V. 

2. The equation above can be understood as defining a matrix model, where the 
general connected matrix-model partition function is given by. 



where M is an n x ?7,-matrix, and ^(M) is the Potential for the model; the 
integral is taken over the space called V above. 

At this point, we note that in spite of all of these generalizations we have done so 
far, we are still talking about a 0-dimensional problem. This is no surprise, for we 
started with a 0-dim problem and changed it in many different ways, but we never 
added any dimension(s) to it. This is the next step, to construct the dimensions. 




(14.2) 




Chapter 



15 



Dimensional Construction Sketch 



At this stage, the question poses that itself is: ''Given the generalization (14.2), how 
does one go about engineering the dimensions in order to extend this problem from 
0-dimensions to d- dimensions?" 

This is a very hard problem and, although it is outside of the scope of this work to 
tackle it, we will give a very rough sketch of what should be done. 

The whole problem of using (14.2) in order to engineer the dimensions from first 
principles is that two limiting processes are involved: the limit where the size of each 
dimension tends to infinite (dubbed "infinite volume limit"), and the limit where the 
number of points in each dimension tends to infinite (called "thermodynamical limit"). 
What we want to do is to "grow" a space whose natural differential operator is given, 
in some sense, by a combination of a derivative(s) and p'($). 

In our case (Chern-Simons), the Schwinger-Dyson equations of motion are symboli- 
cally given by. 



where A i-^ — z 5j, and A, J G so(2, 2). In this case, the operator given by {—i d— AAA) 
is the natural differential operator under study, and it is constructed out of the 
polynomial p'($) = tr A^ = A A A, with the derivative operator being given by 
dA = tr(dA). 

In some sense, the tool that allows us to engineer however many dimensions as 
we want is called D -Module. By definition, a D-module is a module^ over a ring D 
of differential operators. The methods of D-module theory have always been drawn 
from sheaf theory and other techniques with inspiration from the work of Alexander 
Grothendieck in algebraic geometry. The approach is global in character, and differs 
from the functional analysis techniques traditionally used to study differential operators. 
The strongest results are obtained for over-determined systems (holonomic systems), 
and on the characteristic variety cut out by the symbols, in the good case for which it 

module is a generalization of a vector space over a field K when K is replaced by a ring. 



{-i dA 
(-^d(5j) 



AAA) Z[J] = J ; 
5j A 5j) Z[3] = J ; 
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is a Lagrangian submanifold of the cotangent bundle of maximal dimension (involutive 
systems) . 

For more about the use of D-modules as partition functions, we refer the reader to 
[60] (and references therein). In particular, for our problem we want to compute the 
so-called "Airy D-Module" — to learn more about such problems, please refer to [til] 
(and references therein). 

As a final note, we remark on the properties that will be determined by the two 
contours that define (14.1). As we already know, these contours (resp. the boundary 
conditions for the Schwinger-Dyson equations) ultimately define the ranges that the 
parameters (mass, coupling constants, etc) of the theory can attain. More generally, 
when talking about gauge theories, we should realize that, mathematically, the coupling 
constants are the parameters describing the scalar product on the particular Lie algebra 
in question, while, physically, the coupling constants appear in the structure constants 
[defining the Lie algebra under study] and is interpreted as measuring the "intensity" 
of the interaction between fields, [62]. 

This means that different contours will determine equivalence classes of structure 
constants: coupling constants in the same equivalence class determine the same Lie 
algebra, while coupling constants in different equivalence classes will determine distinct 
Lie algebras. 

Therefore, the two contours that yield the Ai and Bi functions and its respective 
extensions (as defined previously) will ultimately determine different Lie algebras, i.e., 
two distinct theories. 

This is not much of a surprise when we think in terms of the results already derived 
in Part II, where the geodesies of the cubic potential (chapter 9.2 of Part II) where 
shown to be significantly different (thus, the group manifolds in question have to be 
just as distinct). 

In particular, in this context of 3-dimensional gravity, one of the solutions is 
particularly interesting, for it vanishes for a certain time, after which it suddenly 
and abruptly starts to grow. It begs the question: does it really mean that, in this 
particular solution, (2+l)-spacetime is non-existent for some time, after which it 
expands quite aggressively? 

These questions of now, still unresolved. 



Part IV. 
Final Conclusions 
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Chapter A. \^ 

Conclusions and Future Work 



To sum up this work, these are a the "milestone results" that should be clear: 

• In general, physically interesting QFTs has multiple solutions, which are associ- 
ated and determined by the contours of integration that render the partition 
function finite, which is absolutely analogous to finding all solutions to the 
Schwinger-Dyson equations of the problem (and its boundary conditions). 

• In order to be able to find all possible solutions, we are forced to extend the field 
configurations to complex-valued fields, or, more generally, to matrix-valued or 
Lie-algebra-valued fields. In turn, this implies that the parameters of the theory 
(mass, coupling constants, etc) are also extended in an analogous fashion. 

• The Lee- Yang zeros of these multiple solutions condense over Stokes' Lines, 
coalescing the solutions and creating "sectors of analyticity" for the partition 
function. Thus, we end up having a "canonical" set of superselection rules in 
order to determine the particular vacuum under study. 

• In this sense, the partition function should be understood as a meromorphic 
function of the parameters of the theory. This is in complete accord with the 
interpretation of the partition function as a particular D-module, constructed 
from the 0-dimensional version of the model at hand, in order to engineer a 
(i-dimensional theory. Moreover, once the partition function is being understood 
as a function of the parameters of theory, which have been "extended" , we can 
focus on its transformation properties, in particular, whether or not it possess 
modular symmetry (SL(2, C)/Z2), which is the automorphism group of the 
complexified parameter space (understood as a Riemann sphere). 

• On a slightly different token, when gauge symmetry is present, the boundary 
conditions of the Schwinger-Dyson equations (resp. contours that render the 
partition function finite) and the gauge symmetry completely determine the 
local system (in the sense of [■!()]) in question and, in this fashion, they determine 
some properties of the Higgs bundle associated with this theory. In particular, it 
is possible to infer the topology of this Higgs bundles, based on the Jacobi metric 
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of the problem. This information is important for the Geometric Langlands 
Duality, where the singularity structure of the gauge connection is of importance. 

• If we are given a pure gauge theory, S = ^ tr(F^), the contours of the partition 
function ultimately determine the coupling constant, which makes its appearance 
in the structure constants determining the Lie algebra of the gauge group at hand 
(analogously, mathematicians understand this as determining the scalar product 
of the Lie algebra), thus measuring the "intensity" of the field interactions. In 
this sense, the structure constants can be dynamically determined by the theory, 
and its different equivalence classes will determine distinct theories. 

Furthermore, in the sense of [57, 63], we can understand the process of starting 
from a "bare" (initial) theory and dynamically determine its parameters (via the 
contours of the partition function or the boundary condition of the Schwinger-Dyson 
equations) giving rise to multiple (final) theories as a quantization scheme. This is 
intimately related to the fact that the quantum corrections, being handles attached to 
the classical solution, deform the classical geometry into a quantum one — a fact which 
is related to the more common Geometric and Deformation Quantization schemes. 

In terms of cobordisms, this can be picture as shown below, where arrows (which 
represent the particular propagator in question) going from the initial theory to each 
one of the final theories represents the quantization of that particular solution: 



Initial Theory 




Topologically Inequivalcnt Final Theories 
Equivalence Classes of points in the Moduli Space 



As for future work, some options are clear, e.g., a better understanding of how 
the boundary conditions of the Schwinger-Dyson equations determine the structure 
constants in a pure gauge theory, and how this affects the different multiple solutions 
of the theory; a deeper realization of which role the Geometric Langlands Duality plays 
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in this framework; and, of course, a more robust comprehension of the 3-dimensional 
gravity problem. On the moUifier front, we hope to be able to tackle more physically 
interesting problems as soon as we have better hardware and software available. 
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Part V. 
Appendices 
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Appendix X A. 

Convolutions and Smoothing 



With the aid of the tools developed below, the smooth approximations of functions can 
be done in a mathematically rigorous fashion. This is useful to justify the statements 
made in this work, where this technology is applied to the generating functional of an 
arbitrary QFT. For proofs of the theorems shown below, see [20] [21]. 

A.l. Mollifiers 

Notation. If U G is open, dU is its boundary and £ > 0, let f/g = {x G 
U I dist(x,(9f/) > e}. Further, let B{0,e) be the ball centered on and with radius e. 

Definition. (Standard mollifier.) 

A mollifier, rj, also called an approximate identity, is a positive ^°°(M"') function. 
The standard mollifier is defined in the following way: 

. Define r] e ^°°(M") to be, 




The constant C > selected so that J^„ r]{x) dx = 1. 
• Ve > 0, set 

Ve{x) = —v{x/e) . 

rj is called the standard mollifier. The functions rj^ G satisfy Vei^) dx = 1 and 
supp(r7e) C B{0,e). 

Definition. (Mollification.) If f : U — >• M locally integrable, define its mollification 
to be, 

Us 3 r = ve*f . 
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A. Convolutions and Smoothing 



That is, \/ X ^ Ue, 

f{x)= / Ve{x - y) f{y)dy = / Ve{y)f{x-y)dy. (A.l) 
A. 2. Properties of Mollifiers 

Theorem A. 2.1. (Properties of mollifiers.) 

1. re^^iU,); 

2. p — > /, almost everywhere, as e ^ 0; 

3. If f E '^{U), then f^ ^ f uniformly on compact subsets of U ; & 
I Ifl^p<ooandfe Ll^{U), then f ^ f m Ll^{U). 



Appendix J—/ 

Entropy Calculations 



In the cases considered in 7.1, the information theoretic entropy can be computed 
analyticaUy. There are 2 types of functions to be considered: 

/.'W = exp{4(^)'}; 

where c is a constant scalar. The first type above deals with the Airy Function action 
for J < (c = ±\/— J ) and the cj)'^ action with < (c = ±a/— yu/2 (7 ), while the 
second type above handles the Airy Function action for J ^ (c = ± J) and the 0^ 
action for ^ (c = ±i ^/jj,/2g). 

The information theoretic entropy has the general form, 



J fi{x) dx + log ( / fi{x) dx 



Therefore, we see that the analytical results are given by: 

(e/2) ^/2^ 



S'(e,c) 



V2tF e + log(v^27r e) 

^/2V (e + c/2e) 
e-'=/2^' e + log(e-'=/2^' e) 
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B. Entropy Calculations 



And, thus, we have closed forms for the entropies of the examples given. The graphs 
below show these entropies: the leftmost one shows S^(e) while the rightmost one 
shows §^(e, c) for a range of c values. 



mrcrnwlbn TlwDralb Emrtfy^ E Ifl 



Appendix 

Parabolic Cylinder Functions 



The differential equation for the Parabohc Cyhnder functions is usually written as, 

-^ + {az^ + bz + c)f = 0; 

where z E C 

The 0- dimensional 0^ action is given by S[(j)] = n4>'^/2 + gcj)^ /A, which yields two 
Schwinger- Dyson equations: 



Perturbative solution: — i 6jZ — J 2, = ; 
Non-perturbative solutions: djZ + JZ — (3Z = Q; 

where (3 = /i/g. The match to the previous form is obtained when a — 0, b = +1, c — 
—f3. Note that there are two non-perturbative solutions, one being solitonic and the 
other one being called simply broken-symmetric. 
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